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INTRODUCTION 

Residual stresses caused by welding result from the nonuniform 
rate of cooling and the restrained thermal contrabMon or non- 
uniform plastic deformation. 

From the zone of extremely high temperature at the weld, heat 
flows into both the adjoining cool body and the surrounding 
atmosphere. The weld metal solidifies under very rapid cooling. 


The plasticity of the hot metal allows adjustment initially, but 
as the structure cools the rigidity of the surrounding cold metal 
inhibits further contraction. The zone is con^ressed and the 
weld is put under tensile stresses of high magnitude Cl) . The 
danger of cracking in these structural elements is great (2) . 

Change in specific volume is caused by the change in temper- 
ature. The following simple analysis helps to visualize the 
action of temperature change in a confined body on stress operation (3) : 

A reduction of the temperature (t) of a confined body by [At) 

^vould decrease the body’s original length ) by C^^ ). To 


restore the original length, 

Hooke's law gives “= 

F. 


. . . 

expansion is c<. ~ . 

X 


a tensile stress (cr~) must be applied. 
, and the expression for temperature 
By substitution therefore. 


At, where >( is the temperature expansion coefficient, 
E 

and E is the body's modulus of elasticity. 


Assuming E = 30 x 10^ psi for steel and‘^= 6.3 x 10”^ per °F, 


the temperature reduction necessary to reach the yield point of 


about 30,000 psi is: 

1 . le^l'l60°F 

E o< 30 X 10° 6.3 



P2 


This indicates a relatively small temperature change in a 
confined body can give rise to an intense stress, and of course 
stresses will result from any temperature gradient regardless 
of whether the body is confined or not. 

In the case of large and complex structural frames, residual 
stresses will definitely be high because of the confinement of 
the individual members. Also the proportions of components have 
significant effect on the stresses . 

As a result of the highly localized heat input, followed by 
the rapid cooling, specific effects are introduced into the 
regions adjacent to the weld. The sudden variation in temp- 
eratures bet^veen the weld zone and adjacent parts, Figs.l § 2, 
can cause non-uniformity in grain structure, and the formation 
of hard Martensite or Bainite, Figs. 3, 4 § 5, which are brittle 
and; with the interplay of thermal stresses, contribute to 
cracking under external factors. The sharpness of the variation 
in the micros true ture can be reduced by preheating the basic 
metal adjacent to the weld, to about 400°F just prior to welding. 
This preheating produces more gradual change in the microstructure 
and thereby eliminate the metallurgical stress -raiser. 

In many projects where stiffness and dimensional stability are 
essential, the removal or at least a significant reduction 
of such stresses is necessary. 
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Oxy -Acetylene Welded Mild Steel Plate (0.2%C) 



Macros true ture - Section Of Weld 5X 

Coarse grains In the heat affected zones (In excess of 1000°C) 
Excessive austenite grain growth near the fusion hoxmdary 



Hlcrostructure 
Parent Material 
Ferrite + Pearl Ite 


500 X 



Mlcrostruc ture 
Heat Affected Zone 
Uouer Balnlte (Cementlte 
Seuareted by wide bands 
of ferrite (E)) 

500 X 


Figure 3 
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Manual Metal ^Axc Welded Mild Steel (0.2%C) 



Kacro structure - Section Of Weld 5 X 

Less extensive heat-affeoted zone corapared to Oiy-Ace tjrlene 
Welding aiaterial subjected to hlsh temperatures for very lonv. 
Small region of coarse grains 


Heat affected Zone 5 OO X 
The large coarsened grains (G) 
probably troostlte or lower 
balnite. are surrounded in 
the boundaries by bands of 
proeuteotoid- ferrite {H), on 
which upper bainite (J) has 
nucleated. 



Figure 4 
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Some o£ the patterns of residual stresses are understandable, 
especially when a metallic solid is exposed to a homogeneous 
treatment . Various theories that depend upon layer removal 
(4-6) have been developed to analyse uniform stress distributions. 

The decision was made to release the embeded residual stresses 
by removing successive layers off the metallic surface using 
diluted nitric acid solution or caustic soda solution. Although 
such a technique is time consuming, yet it is preferred to the 
mechanical methods where excess heat and vibrations can affect 
the pattern of the existing residml stresses. 

Fusion welding is generally used because it is rigid and simple, 
as compared to rivetting or bolting techniques. In fusion 
welding tensile stresses of high magnitude are ali^ays developed 
within or near a weld(7) and the danger of cracking in these 
structural elements is great. 

Thermal treatment is currently applied to relieve such stresses. 
The technique of inducing mechanical vibrations, through a 
variable speed eccentric motor, has been tried with simple 
metallic shapes to reduce or eliminate residual stresses developed 
by welding (8-10) . The effect of such a technique on complex 
structures, however, is not clear and needs investigation. 

One major project affected by fiasion welding is the Laivrence 
Livermore Laboratory's Shiva project in Califomia(ll) . 
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Stablein, Davideiikow and Slievadin., Letner,^^^^ and Rrisch 
and Thomsen, all tried to deduce mathematical equations for pre- 
dicting residual stresses in bar specimens, applying the deflection 
technique. The most suitable analysis is the one developed by Frisch 
and Thomsen. If a straight metallic bar of rectangular cross section 
is siibjected to such loads that stresses are induced in one of its 
surface, the stressed surface of the unrestrained bar will become curved 
upon release of the load. Furtliermore , if the induced stresses are 
tension stresses and are constant at any particular depth in the stress- 
ed layer, the bar will assuiae a concave curvature of radius R, 

The analysis deals with longitudinal stresses only, and it has 
been assumed tliat the effects of possible transverse stresses, if 

present, are negligible. ORIGINAL PAGE IS 
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Vdiere: 

Ojj = Original stress in the infinitesimal layer A hj^ which is n 
arbitrary infinitesimal layers of thickness A h below the 
surface of the original bar 
E = Modulus of elasticity 

2L = Central portion of the bar between supports 
H = Initial depth or thickness of test specimen 
Ah 2 = First layer to be removed 

A £ = The change in deflection after the first layer is removed 


Successive layers of thickness Ah are to be removed from the bar 
until hg is completely removed. From the observed changes in deflec- 
tion of, the bar, it is possible to reconstruct the original stress dis- 



The above relationship can be simplified to give 
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K = A E = A constant for anyone material 

d£ = The slope of deflection as a function of the depth of the 
dh 

layer ' . _ . 

Also it may be \\TriT:ten in a logaritimuc' form while replacing a by 
n, the true response of the residual stress (psi) 


n = Bq + B;^Xjl ^ B2X2 + B3X3 3 

or ■ . , 

Y = bo + biXi + b 2 X£ + b 3 X 3 + e 4 

\ilhere: 


Bq, Bj^, B 2 , i B 3 = Parameters 
bo, bp, b 2 , § b 3 = Estimates of the parameters 
e - Ej<perimental error 

Xp, X 2 , and X 3 = Logarithmic transformation of }\2 and df 

dh' 

Y = The observed residual stress on the logaritliniic 
scale 

The stress in the first layer can be obtained by obtaining the 
slope of the deflection versus depth- of- layer curve at the surface and 
letting A h 0 , Thus 


Ol -B.hZ , df , 

3 ^ dh -^surface 


5 


IVhere: 


dh 


surface = Slope of deflection as a function of depth of 
layer at h=0, or surface of bar 


Revealing the pattern of the residual stress distribution for each 
of the geometrical shapes under consideration will help to develop a 

new set of analytical formulas to represent such stress patterns as 
fmmd Tn romnl PV f-raine structures bv welding. 



Original page 
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(Hi-2h +h +Ah.h/2 




Figure 7 


Assume the hody of a "bar is formed of infinitesimal layers 
of thickness Ah each. 

Let he the stress generated hy welding in a layer (n) 

below the original surface of the bar and of thickness Ah, 

n 

Following the removal of all layers above 


original stress 


in that layer will be changed to 


This change is caused by the redistribution of the stresses 
following the removal of (n-1) layers above the layer (n) and 
in order to restor^the balance of stresses across the bar section 
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In order to visualize the effect of the removal of any 


specific layer, at depth h^, upon the stress cr^ in the 

layer « “ 

Let o-;; be the stress in the layer Ah at the time this 
m m ' 

layer is removed. 


The removal of the layer ^h causes two changes upon the 

m 

stress cr~;j“ 
n 

1- Inducing a stress CT~^ because of the removal of the axial 


force in the layer J^h 


m 


2- Inducing a stress' CT^g because of the removal of the bending 


moment effect of the force removed in the layer Ah 


m 


Therefore, the original stress cT*^ in the layer A h^ 

caused by the welding is composed of i- cr^» which is the 

n-1 

stress in that layer at the time of its removal + (o — + O 

Q nA nB 

which is the total stress induced in that layer because 

of the removal of all previous stressed layers. 


t « 




n-1 

+ r< 


0 



( 6 ] 


If Pjjj = The axial force in the layer removed, it balances 
the remaining axial force F 

O 

P = - P 
s m 

Assuming the width of the bar = b 

^m ~ ^ cross sectional area of the layer 

~ Ah^.b (73 
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The remaining axial force F_ is worked from the balancing 

s 

axial stress cJ“. appearing at layer n. 

iU% 

"s = - \)-'= 

. „ Ps I'm 

‘ ' 'gfi. = - 

— ^ (H-h^).b (H-h„).b ’ H-h^ 

The second reaction is the bending reaction 
M„.C„ 

c— ■ = (10) - 

nB j 

■^s 

I-There j - 

I ^ — Clla) 

® 12 

= The distance of the layer n from the HoA. of the 
remaining section 

^ H-2h^^-h„-PAh„ 

2 

M„ = The bending moment in section remaining 

5 

= P . Its arm from the DJ.A. 

( 12 , 

“■ 2 

Substituting (lla) , (11b) and (12) in (10) 

-P . . 

' S 

b(H-h 
^ m 

12 


0 
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And substituting for from equation O') 

• _ (;i3b) 

=S~ (H-hj,)3 

Stress at the outer fiber 

CT-* == p. ^ (14) 

li'Hiere 

E = Modulus of Elasticity 

" Strain measured using electric strain gages 
The total change at layer n, due to the, removal of ^ 

trn^ + 'HIb 

From equations ' (9) ard (13b) 

_ m* 3 <7~^ ■ ^ ^ H-2hy^+h^-t 

“1 H - h„ (H - hjj)3 


9s. OV 


^ nl 

^ 3A\(H-h„+Ah^)(H-2h„H-h^+ah„ 
H - IV (H - 


(16b) 


And for the outer fiber 




And A h^ = 0 


C7-^. Ah^ 

3 in* Ah^(H-hj^+ Ahjjj) (-H+hjj^) 

H - h 

m 

(H - h^)3 


= _ ^ 



H - 

(H - lv)2 


(17a) 
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c— 6 

~ ~ 


3Ah^^(H-hjjitAhjji) 



(H-h^)2 



E. 



3 &h„(H-h„+ Ah^) 



(H-H^)2 

Prom equations (16b) § (17b) 

. ( o — = _ 


X 

Ah 

“m 


H- 

-^m 




3Alv(H-h^+Ahj,)(H-2h„+h„+Ah^) -| 


(H - J 


PIO 


(17b) 


(17c) 


- 3Ah„(H-h„+/ih„) 


(H - 




O’ 


(H~ly) + 3(H-h^+Ah^)(H-2h^-Hl^4-£>h^) 

(H-hj,)^ - 3(H-h^+Ahj,)(H-h„) 


And 


nl 


E. €. 


O’ 




(18b) 


This equation can be put in the form 



Inhere 
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~ ^ ~~ ^tn 

■B = H- h^ + Ahjj = A + Ah„ 

C = H - 2h„ + h„ + Ah„ = A t 2(\-h„) +Ah^ 


Since we need to Imow the change in stress at each level 
for each layer removed, A'^q s.nd. are introduced as the 

change in strain at the outermost fiber and the change in the 
stress at level n respectively 


From equation 


(18b) 


y 

I ^ 



= E. AC^ 


+ 3tH-h^+Ah„)(H-2h„+Ah„+h„) 
(H-h^)2 - 3(H-h^+Ah^)(H-h„) 


CIS) 


Procedure »- 

1- Remove a layer by etching 

2- Calculate the stress in that layer using equation (17c) ; 

It will be the residual stress in that first layer 

3“ I'/hen removing each subsequent layer, the stress in each 
layer is the sum of the residual stress in the layer plus 
the change in the stress due to the removal of the previous 
layers . 

Equation to be used is equivalent to equation (6) 


i.e, , cr~^ ~ 


— 

nR + ^ 'Acr^i 
1 



n-l 

2T A 


(20a) 


Or 


1 


nl 


(20b) 
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This means that equation (19) is used to find the changes 
in stress at any level. Sum them. And upon the removal 
of layer n, equation (20b] is used to find the residual stress 
in the nth, layer. 

Example 


H = 0.2505 inch 

First layer removed Ahjj^ = O.OOO6 h^^^ = O.OOO6 inch 




= OoOOOO 


4 X lO"^ in/ in 


E = 29 X 10 psi 


Prom equation ( 12 c), 

29 X 10° X 4 X 10“° 


Cf — = - 
m 


0.0006 


(0. 2505-0. 0006) 


3 X 0.0006(0.2505 - 0.0006 + 0.0006) 

(O.2505 - 0 . 0006)2 


116 X 0.062 
0.0005 


+ 24,000 psi 


Since this is the first layer, this is the residual stress 
in that layer. 

For the change in stress in the next layer 
Ah^ - O0OO59 


hj^ = 0.0059 + 0.0006 = 0.0065 

Prom equation (l 4 ) 

A = 29 X 10^ X 4 X 10~^ 


0.2499^ + 3 X 0.2505 X 0,2440 
0.2499^ - 3 X 0.2505 X 0.2499 



?13 
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- 29 X .14' X 


0«062 + 0.183 
0.062 - 0.188 


= - 226 psi 


Second layer removed i- 


= 0.0059 , iTji = 0.0065 , Cq = 1 X 10' 


prom equation (12c) 


29 X 10^ X IX 10"^ 

0.0059 3 2 c 0.0059 X 0.2^99 

0.2440 (0.2440)^ 


1.74 

0.003 


= 580 


And from equation (20b) 

-(■- - 226 ) 
r= 8 O 6 psi 

which will be the residual stress in this second (nth) layer. 



The basis for the calculations of welding heat originated in the work of 

) 

Ibberts And Rosenthal [16,] 7) .Various assumptions were necessary in the derivation 

of the mathematical formulae, including the possible radiation and 

convention losses (18). Heat is lost by conduction into the adjacent 

parts of the metallic bars; in addition, the bars lose heat from its 

surfaces by convection and radiation, which was proved previously to 

be negligible (19) . Radiation loss can be calculated from "Stefans 

Law" to be about 1.8 cal/square centemeter or 7.5 watts/square centemeter, 

if the temperature is 800°C. If the remainder of the surface reaches 

200°C, the radiation loss will be 0.3 w/sq.cm. Convection loss will 

be of the same order (20). 

In a theoretical analysis it is necessary to know accurately the 
amount of heat reaching the work-piece. Efficiency with Vee-Groove 
and a good arc may reach 90%, while it may drop to 30% with simple 
melt-run on a flat plate with poor flames. 

In their mathematical analyses the above authors aimed at throwing 
some light on the effect of the material characteristics and welding 
conditions on the temperature distribution, but it is not intended to 
supplant experimental measurements, they related lieat input to vd/4 
IVhere V = welding velocity 
d = fused width 

theoretical diffusivity of the parent material 
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Example: 



oC e.g.s units 

arc -VO It 

thermal 



conductivity 

Mild steel 

0.082 

30.4 

^ 0.10 

Aluminum Alloy 

0.88 

25 

0.45 


If vd/4 is ^1, the heating is efficient with not much losses 
Also Q = 81cTj^Cl/5 + 4d/4 ) 

I'/hexe Q = rate of heat input to plate Ccal/sec/cm) 
h = thermal conductivity 



PROJECT DESIGN 



PROJECT DESIGN 


The suggested investigation i'?ill 'cover structural steels. Type 
ASTM A36, \d.th flat' and V/- cross sectional shape, to 1/4” 'thick. 
Reinforced box geometric structures will be built using "the semi 
automatic solid v;ire arc welding technique, to confonn to AV/S Bl.O 
standards, ( 21 ) using Hobart wire #28, AWS E70s“6, size 0.035 as filler 
metal (Consumable wire AWS A5-18 for mild steel electrodes) with 
maximum current of 100/225 an^jS and msbdmum voltage- of 19/28 VQl'ts„ 
Edges will be brush free of mill scale, rust and other foreign matter. 
The base metal will be dry and free of all oils, paints and foreign 
matter. Temperatures will be checked by -the use of temperature indi- 
cating crayons . 

All joints will be of the Groove-fillet type. Tlie filler metal 
will eliminate possible variation in the coefficient of linear expan- 
sion of the parts and prevent complications in "the stress analysis. 

The geometry of the -structure will be studied closely . Marks 
•will be made at certain distances around 'the ^'^elded zones and accurate 
measurements will be taken. 
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Electric strain gages will be fixed on the metallic surfaces 3 protect- 
ed by suitable adhesives against the acid- treatment which is to be used 
to remove successive layers for releasing residual stresses. 

Several frames of the same structural shape are to be prepared 
and used for the various phases of the study, that will include: 

Ca) Coirrparison of the original and welded shape and dimensions 

f 

(b) Sectioning and releasing of the residual stresses by layer 

removal with diluted nitric acid solution, a technique that 
has been developed by Waisman and Phillips. ( 5 ) This layer 
removal method is most useful in determining the residiml 
stresses in flat plates .or beams in which the stress varies 
with the thickness, particularly if the stress gradient is 
sharp. This method consists of rem.oving thin layers in the 
region of interest and observing the dimensional chages in 
the part from xdiich the layers were removed. The average 
stress in each layer can be calculated and a plot of the 
residual stress distribution obtained. The same chemical 
etching technique has proved suitable and effective for other 
shapes including the cylindrical, solid and hollow geometries .(22 

(c) Exposure to vibrational conditioning as a stress relief 

' technique, by attaching a stress relief vibrator to a central 

area of the frame, in a working cycle of fifteen minutes. 

Dasgupta, Baker and Weiss (2^ have studied the vibrational technique 
to relieve residual stresses. They learned that dynamic surface stresses 
are induced in or near a weld. 

The basic stress relieving method enployed will be that is referr- 
ed to by The Stress Relief Engineering Co. of Costa Mesa, California 
and The Rigging International Of Oakland, California, as The Formula 
62 System. 
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The procedure recommended is to relieve stress immediately after the 
frame is completely v;elded by placing it on level rubber pads and run- 
ning vibration tliru its 15-minute automatic mode. Tlie technique is a 
non- thermal method of relieving the residual stresses in metals or metal 
weldments without distorting or dianging tlie tensile strength, fatigue, 
or yeild poiiit of the work piece. It is necessary to prove that tlie 
vibrational technique is reproducible and capable of being monitored 
and effective at all critical locations witliin a part being so processed. 

In the thermal stress relieving process, energy in the form of heat 
is introduced into the worlcpiece to reduce stresses. In tlie mechanical 
process energy is introduced into the ivorkpiece by means- o£ vibrations 
in a low frequency range of 0 to 100 q3s, which will generate enough 
dynamic stresses to interact with the high intensity surface tensile 
residual stresses, and at the same time eliminate the side effect of 
vibrations using up a large fraction of the fatigue life of tlie'p'art 
being processed. It is important for this investigation to establish 
relationships among stresses, forces, amplitudes and accelerations de- 
veloped during vibration. 

Relationships bet^^^een force and. deflection characteristics should - 
also be experimentally investigated. 

An accelerometer can be used to monitor acceleration of a vibrating 
part. Electric strain gages can be used to measure about half the yield 
strength of the material as an indication of the level of vibrations 
to be applied for stress relief. 

The resistance strain gage can. be used in conjunction -i^rith -the 
accelerometers to characterize the vibrational act for complex shapes 
and to monitor the vibrational stress relief operation. It may also 
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be liclpful in indicating vibrational stress relief v;hen oriented parrallcl 
to the direction of maximum residual stresses. 

Combining measurement and analysis to characterize the vibrational - 
behavior of a given shape, and monitoring the vibrational stress relief 
treatment, reCiUires a spectrum analyser. Such an instrument exists at 
N; C. Agricultural and Technical University, similar to the more expen- 
sive Fourier Analyzer, except .that it needs continuous vibrational 
source applied to the structure. The accelerometers and strain gages 
are affixed to the part, which is then subjected to a continuous vibra- 
tion. Signals are processed by the analyser and the modes associated 
with the natural frequency, characteristic of the part, are determined 
and displayed in animated isometric form on a cathode ray tube and 
photographed. 

The frames, vibrationally conditioned, will be cliecked for stress 
relief in tivo ways: 

Caj Tolerances, squareness and parallelism will be measured and 
compared with a similar frame untreated by vibrations, 

(b) A full stress investigation by sectioning and layer removal 
will be made as explained before. 


PART ONE 
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Iffork on the project started in April, 19?8 and has continued 
for six months.. The first two and one half months I'rere spent 
purchasing .equipment , strain gages,' accessories and steel 
■bars. A stand for etching the metallic bars and a tray for 
holding the acid solution were d.esigned and built of plexyglass’ 
sheets, 3/8 inch thickness, cemented together by ethelyne 
dichloride, Fig, 8 

Carolina Steel Co. of Greensboro, W.C, through Mr. Robert 
Sherman, Vice President for Engineering has encouraged the 
project, showing support by offering to cut, weld and recut 
the steel and boxes free of charge when I purchased the required 
steel from their company. 

Electric strain gages were supplied by the BLH Co. They 
are designated SR-4, FAE-50-35 SX (350 Ohm) with a length of 

0.500 inch (Gage Factor 2,07) and 0,250 inch (2,0^), The gages 
are encapsulated with a 0,001 inch polyimide film that nrovides 
protection for the sensing element during insulation handling. 

The film also promotes better long term stability as the foil 
grid is protected from air borne contaminants or fingerprints. 

The gages are self temperature compensated j They automatically 

cancel unwanted changes in gage resists,nce, due to temperature, 
so that only stress produced, strains are measured. The (^gree 
and accuracy of self temperature compensation is defined in 
terms of a chord slope temperature coefficient; apparent strain 
is (tJ in/in) per degree F temperature change over a specific 
temperature range on a material having a specific coefficient of 
expansion over the same range. 
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Recommended temperature range of these gages is -100 to 
400 °F. The gages are made of constantan on polylmide carrier. 
They were provided with nickel-clad copper-ribbon leads soldered 
to the gage tabs (one inch long) with 425 °F tin-silver solder. 

To prepare the metallic surface at the spots where the gages 
are to be fixed, the steel is smoothed with a sanding disc 
fixed to a portable rotor. . The surface is cleaned with Methyls 
Ethyl Keton. The surface is next roughened with sand paper and 
re-cleaned with MEK. 

• A layer of EPY-500 adhesive- is used to adhere the gages to the 
metallic surface » The adhesive must be cured at about 1?0 
for 24 hours using electric heaters under slight pressure of 
about 20 psi, using pads and wooden blocks wrapped with string. 
This type of adhesive can stand a strain of about 5 % at room 
temperature. It has an excellent moisture and chemical resistance 
and is useful over a temperq.ture range of -450 to 600 °F. I'fhen 
the wrapping is removed, the gages are checked for adhesion, 

A gage tester, is used to check on gage adhesion and insulation'. 
Readings should show less than 1/10 % deviation from a resistance 
of 350 Ohm, which shows that the gage is within the range that 
the bridge can balance. Lead wires with polyvinal chloride 
tephlone are next soldered in place, Mylon plastic tape is used 
to isolate the wire ends off the metallic surface. 
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A barrier D which is -a two component -semi-rigid epoxy system’, 
is spread on the gages to protect them against water and chemicals. 
It has a temperature range of -100 to 225°F. It has. to be cured 
at a temperature of about 150°P for one hour using the electric 
heaters. The layer is thin in order to not stiffen the gage. 

The recording equipment is composed of a portable BLH digital 
strain indicator, model 1200A solid state, and a switching and 
balancing unit model 1225 that -can provide up to ten input channels. 
The recorder has an accurate setting of gage factror of 0.001 ’ 

and a four-way display of +10,0.00 inch per inch, 

•Before etching takes place, '-the gages are- connected to the 
s.train indicator through the switching and balancing unit. The 
gage factor is adjusted and the indicator is balanced for each 
gage. To compensate for the temperature effect, three leading 
wires are soldered; two to one side and the third to the other 
side of the gage leads, Fig.- 9 

Structural A3^ steel bars, hot rolled, 6 inch wide by i inch 
thick are sawed, to a length of 3 feet each. All cutting is 
done by means of a saw. No flame cutting is used. This 
eliminates any heating other than the welding itself.. 

One bar is tested as received .in order to investigate any 
possible residual stresses .while forming. The analysis showed 
negligible residual stresses in that bar. 

Square boxes, 3 feet by 3 feet, Pig.io were formed by welding 
the- hot rolled bars together either by continuous real 
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intershield welding technique using 26 volts, 275' amps., 3/32 
inch hare wire and powder flux, or by an AC i^lder using 225 amp 
5/32 inch rod.- Such shapes simulate some of the conditions 
in built-up structures. 

The sides of the box used are sufficiently long to obtain a 
unifrom state of stress distribution. Strain gages are fixed 
along the bar and simultaneous measurements are taken after each 
etch. One gage is fixed across each bar. 

Future studies will cover thicknesses up to one inch and 
the width of plates up to 15 inch. Such study will reveal any 
plate distortion that may take place beside the actual residual 
patterns generated in such plates. 

At present a set of tests is being performed on ^ inch bars 
welded together in various ways to form square boxes and then 
sawed and analyzed as has been done before, 

A third group of tests is planned to weld unequal sections 
together in order to study the effect of various widths when 
.welded together and if such a structure may develop concentrated 
residual stresses. 

Sectioning and releasing residual stresses by means of layer 
removal, was performed using nitric acid solution of 30 ^ 
concentration. This technique has been .developed and proved to 
be most useful in determining the residual stresses in flat 
plates as well as many other sectional forms 
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Etching precautions include wearing respirators for toxic 
fvunes and gas, rubber gloves and rubberized apron on top of 
lab coat. The stand and the acid tray, are supported on a table 
under a hood with an outlet fan running continuously. 

Nylon scrub brushes and paint brushes are used to continuously 
brush the metallic surface with the acid solution; because 
of the inclination of the supported bar the residue drops into 
the acid tray. 

The bar is etched, After 24 hours measurements of strains and 
thicknesses are taken, to be followed by a second etching and. so on 

Magnetic surface temperature thermometers are fixed to the 
back side of the bar while etching. No excess heat is registered 
{170 to 180°P). 

3 inch and 6 inch deep throat micrometers are used to measure 
the bar thickness at various locations. 

The average stress in each layer is calculated using the 
theory explained and presented in this report, together with a 
developed computer program, a copy of which is enclosed. 

Results of the residual stress distributions in the bars 
cut from the welded straight and. box shapes are shown in Pigs. 

11 to 18 , High tensile residual stresses are found at the weld 
(about 25,000 psi) dropping sharply to compressive stresses in the 
range of 3,000 to 5,000 psi which spreads for some depth through 
the bar section. Surface stresses change gradually along each 
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bar, from maxiimni tensile of about 25,000 psi to maximum 
compressive of about the same value at a distance of about 7 to 
8 inches away from the weld. 

The similarity between the patterns. of stresses obtained 
from the flat welded bars and the box welding indicates 
that the _-^ress distribution in various shapes can be obtained 
approximately by superposition. These stresses can, therefore, 
be considered as the cumulative sum of the thermal stresses set 
up by welding. 

In order to judge the reaction on the strain gage fixed across 
the bar, other tests are planned with various bar widths and ' 

■with more, than one gage fixed across the span. 

Both tests No. 1 & 2 are performed on. bars sawed from a set of 
three bars, 3 feet each, welded together end to end as shot>m'V'! 
in Fig. 10. Nelding is carried out with the #3 AC Jet Welder 
mentioned before. 

Test No. 3 is performed on a bar sawed from a square box welded 
together using the Inter Shield Technique, 

Test No, 4 is performed on a bar sawed from a square box welded 
together using #3 AC Jet welder, 

With these results and with what may be revealed in the further 
built-up shapes planned and ■under study at present it will be 
possible to estimate the general stress distribution in welded 
built “Up shap e s . 

This is essential before going through the second phase of 
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of the study that will deal with shaking the various shapes of 
built-up shapes, registering the spectrum and analyzing it through 
a Fourier Analyzer using both a load cell with strain gages, and 
an accelerometer with a suitable amplifier. It is possible to 
register such spectrum and then feed it into an analyzer using a 
suitable computer program. The structure is to be shaken for about 
15 minutes at the freq.uency that shows maximum amplitude, and 
then will be tested by a complete residual stress analysis to 
investigate how the welding stresses can be relaxed by the 
vibrations. 

Standard tensile test specimens ^ inch x i inch are prepared 
lengthwise and crosswise from one of the bars as received. Strain 
gages are cemented on the specimens and calibrated against load. 
Results are as shown in Pig. 19 The calibration was carried out 
on a 10 ton tensile testing machine. 

Fatigue testin'g specimens are prepared from one side of a 
xvelded box as' near to the weld as possible. These specimens 
are being tested at the present time. Similar analysis will be 
carMed out on specimens from a box after vibration, in order 
to investigate how such vibrations may affect the fatigue life 
of the steel. This is considered as an essential study especially 
where the steel is to be exposed to dynamic loading in practice. 

Definite conclusions cannot be made at this early stage of 


the study 
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The objective of the first phase of the pro;?ram“ is to determine 
experimentally the magnitude and distribution of stresses 
lengthwise and across the thickness of each member. The use 
of electric strain gages eliminates the need for straightening 
steel merabers following the welding. 

1»3 by volume diluted nitric- acid solution is used to etch 
successive layers off surfaces opposite the gages, The rate 
of metal removal is 2 to 3 mills/l5 minutes. Dimensions and 
strains are measured after a period of 24 hours following each 
etching process. This arrangement permits measurement imder 
stabilized conditions. 

The stresses are next calculated using the theoretical 
approach that we developed 

6 inch and 4 inch wide by % inch thickness hot rolled steel 
bars are used in the first and second phases of the study. 

The bars are cut to a length of 3 feet by means of a saw. No 
flame cutting is used, eliminating any heating other than that 
from the welding itself.. 

Bars are also tested as received to investigate any possible 
residual stresses from the milling. The sbudy showed negligible 
stresses In such bars. 

Square boxes 3 ft x 3 ft are formed by welding which simulates 
some of the conditions in built-up structures. 

Future investigations will study multi-pass welding, weld 
widths up to 15 inches, and unequal section welds to determine 
the effect of such connections on concentrated residual stresses. 
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Current results indicate high tensile stresses at the’ weld 
of about 25,000 psi, dropping sharply across the bar' thickness 
to compressive stresses of about 3^000 to 5,000 psi. Surface 
stresses shovred more gradual changes along each bar, from the 
maximum tensile of 25,000 psi to maximum compression of about the 
same magnitude at a distance of 7 to 8 inches away from the weld. 
Some flat bars were welded end to end; a length of 3 ft. 
sax-jed with the weld central. Patterns of stresses obtained 
from this study are similar to those obtained with box welding. 
This similarity may indicate that the stress distribution in 
various shapes can be obtained by superposition. 

How to relieve such stresses? Traditionally, residual 
stresses are relieved through thermal treatment, suitable for 
simple and symmetric structures. Complex and heavy structures, 
however, cannot be thermally stress-relieved. 

The technique of inducing mechanical vibrations in the x<ireld 
to reduce or eliminate residual stresses has been tried with 
various metallic shapes. Consideration must be taken of the fact 
that a continuous application of vibrations or vibrating 
near resonance, may lead to metal fatigue, and a possible 
spreading of any erabeded cracks. For this reason it is 
planned to vibrate ‘at a frequency within 100 cps, which would 
generate enough dynamic stresses to interact with the high 
intensity tensile stresses and at the same time eliminate 
metal fatigue. 
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In the second phase of the study which is in progress at 
the present time, a load cell is fixed to steel structures 
similar to those already analyzed for stress patterns, A 
hammer is used to strike the cell generating an impulse to the 
system, Pigs. 20, ?-i . It is important to observe the frequency 
response of the force input before proceeding with the test. 

The signal will be treated in either one or both the following 
two ways;- 

1“ The signal will be sent into a channel of a Fourier - 
analyzer. The acceleration response is measured through an 
accelerometer, amplified and sent to a second channel of the 
analyzer. Figs, 22,23* 

The impulse should have a low enough freq.uency input to excite 
the system sufficiently at low frequencies, A softer or 
harder hammer may be indicated by the observed frequency content 
of the impulse. 

2~ A second method, which is being applied at the present 
time is to read the data of the pulse and feed this information 
into a Dec, 10 computer available on campus, through a suitably 
developed program that uses data digitized manually and stored ■ 
on a magnetic disk. It is designed for the joint analysis 
of two random, stationary, ergodic processes containing 
significant noise. The frequency analysis consists of Fourier 
transformed probability functions and statistical smoothing. 

The fundamental utility of the program is in the analysis 
of frequency response characteristics. 
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The routines used include* Power spectrum density, cross- 
spectral density, Transfer function analysis, Coherence 
analysis, Hanning window smoothing and Parzen lag smoothing 30t3l 

The structure will be shaken for about 15 minutes, using 
a mechanical shaker, Fig. 2^ at the frequency that will show 
maximum amplitude, according to the previous shock study. 

The third phase of the program .-will be to test the vibrated 
structure and carry out a complete residual stress analysis, 
as in phase one, to reveal the effectiveness of the vibration 
technique on the embeded residual stresses. 

It is also important and necessary to prove that the vib- 
rational technique is reproducible and capable of being monitored 
and effective at all critical locations within a part being so 
treated. 

Fatigue test specimens will be prepared from the welded 
structures and tested for fatigue life, and similar analyses 
will be carried out on specimens to be prepared from similar 
structures after vibration. 
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Future Work Plan 

Study multi-pass welding with har thickness up to 1 inch, 
weld widths up to 15 Inches, and unequal section x^elds to det- 
ermine the effect of such connections on concentrated residual 
stresses, 

A load cell x-fill be fixed to steel structures similar to 
those already analyzed for stress patterns, A hammer will 
strike the cell generating an impulse to the system. The signal 

will be treated in either one or both of two waysj- 
1- Amplify load cell and accelerometer signals and send to 
txfl-o channels of a Fourier analyzer, A softer or harder hammer 
may be indicated by the observed frequency content of the impulse. 
2- Read the data of the pulse and feed the information into 
the DeclO computer through a suitably developed program, 

The structure will be shaken for about 15 minutes, using 
a mechanical shaker, at the frequency that will shoxv maximum 
amplitude, according to the shock study. 

Next the vibrated structure will be analyzed for stresses, 
as In phase one, to reveal the effectiveness of the vibration 
technique on the embeded residual stresses. 

Fatigue test specimens will be prepared from the welded 
structures and tested for fatigue life, and similar analysis 
will be carried out on specimens to be prepared from similar 
structures after vibration. 
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A second set of steel boxes was formed by welding bars 

f 

^ jc ^ inch section. Details of shapes and dimensions of 
these boxes, referred to as tests Ho. 5, 6, and 7, are 
given in Pig. 25 The results of the stress distribution 
measured as before are presented in 'Figures 26 to 31. 

Test 5 was stopped short due to an apparatus failure, 
yet the general trend of its results together with those 
of test 7 are similar to the previous \results of tests 1 to 
4 of the first report. 

The second phase of the study developed a set of computer 
programs to analyse, the load pulse and the acceleration 
response that were obtained when steel frames similar to 
those studied for stress analysis were secured to the ground 
through brackets, A copy of these programs, together 
with the theoretical analysis and a key to applying these 
programs are presented later in this report. 

Vibrational characteristics of a part have to be analysed 

experimentally. If the part Is vibrated near its natural 

frequency, vibration amplitudes become quite large, even 

though the input forces or displacements are relatively 

small. Correspondingly large dynamic stresses can be 

developed, which interact with residual stresses, A typical 

set-up for impulse testing the steel frames is shoim in 
Pig. 32-1 
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A load cell provided with an electric strain gage, and 
an accelerometer were attached at opposite ends of the top 
side of the frame. A plastic hammer applied a force to the 
load cell at point X, Pig,i2-I and the signal was sent into 

e 

channel 1 of an oscilloscope. The acceleration response 
was measured at point Y using the accelerometer and sent to ■ 
channel 2 of the oscilloscope, A photograph was taken of 
the display. 

Results of treating a frame similar to the one tested 
for stress distribution ITo.? are shown in Pig. 33 a blow- 
up of the photograph was prepared, ■ from which 

measurements of both the load and acceleration frequencies 
were taken. These values represent the X-data and Y-data 

«r 

for the computer programs. The output of these programs 
give statistical analysis of the frequencies, while the 
Trans, results give the frequencies in cps against the 
magnitude of deflection. Results of the Trans, values are 
presented in Pig , 34 showing the prominent periodic frequency 
component that may help reduce the stresses and may reveal 
the frequency of the structure. 

The frequency chosen for vibrating the frame in testS 
was about 1.62 cps. An K3 model CIO vibra.tion exciter, 
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tilted to a horizontal position was used. The frame xvas 
connected to the exciter with a. solid rod, as shotm in 
Pis* 32-11 ■ > This arrangement was necessary because the 

heavy weight of the steel frames prevents placing them on 
top of the vertically positioned exciter. The exciter used 
has a frequency range of 5 to 3000 cps that can be controlled 
through an attached panel. The vibration is accomplished 
principally by means of an induced electromagnetic force 
which is localized in the region of the basic structure 
or body casting. 

The frame imder study was vibrated at l62 cps for 15 
minutes. It was removed, the sides were sawed apart and 
a stress analysis was performed on the too side that was 
exposed to the heaviest vibration. 

Results of this study are presented under testS, in 
Pigs. 36 § 37. , These results indicate that the vibration 
affected the stress distribution at the critical location, 
i.e., at the weld. Vibrating the box developed sufficient 
dynamic stresses to reduce the intensity and general distri- 
bution characteristics of the residual stresses generated 

» 

by the original welding. 
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Quantitative data from measuring, residual stress -Dattems 
■both after welding and after vibration are needed in order 
to state some definite conclusions. Also it is iranortant 
to determine how vibration treatments effect the various 
steel structures to relieve welding stresses. 

Although experiments performed up to the present time are 
limited, the results demonstrate that applying a suitable^ 
vabration amplitude to low carbon steel weldments signi- 
ficantly reduces peak residual stresses in the weldment. 

Such vibration seemed to smooth out the pattern o'f residual 
stresses by reducing the localized high stresses due to 
welding. 

Such findings were expected due to the fact that residual 
stresses are the result of Internal nonuniform ‘Dlastic defo- 
rmation and that the vibration reduces the disorder by 
changing the dimension of the structure. The procedure, 
however, is expected to be more complicated with complex 
shapes and multiple modes of vibration may be needed to allow 
all junctions of the structure to be exposed to the suitable 
vibrational amplitude. 

The effects of the vibrational technique and the corres- 
ponding induced cyclic stresses on the fatigue of steel is 
still under investigation and no conclusions can yet be 
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presented at this stage. 

The boxes formed and studied are of comparatively simple 
shape. The following studies will cover more coraolex shares; 
pulses will be induced at various locations and analysed. It 
is hoped that through such studies, a standard type of 
treatment can be formulated. 
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The ■welded steel boxes studied during the period from April 15 
to November I5t 1979 "were formed of bars 4, 6 and 8 inches wide 
arid and % inch thick, They had various shapes. Table 1 , and 
other complex shapes are planned for the coming study period, 

Figs,38»58 show the results of the tests conducted on seven 
sets of boxes. Each set of boxes was sub;Jected to three phases of 
testing. The first phase of the study examined the residual stress 
distribution across the thickness of each bar and at various locations 
along the bar directly after welding. Both the theory and the tech- 
nique for such analysis together with the developed computer program 
for stress calculation have been presented , - 

The second phase showed the results of impact testing using a 
load cell, an accelerometer and an oscilloscope. The study was 
performed on a second box identical to the one examined in the first 
phase. This phase also analyzed the pulses, applying the computer 
programs oresented at the 'end of this report & the chart of 
displacement vs, the frequency of vibration, .,A specific frequency 

was chosen to be x«rithin 100 cos that showed a peak displacement. 
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For the third t>hase of the study the second box was fixed on the 
floor facing the mechanical vibrator, A steel rod of L-section, 
with welded studs at both ends, was used to connect the box' to the 
vibrator. Other designs will also be used to suit other shapes of 
boxes for future investigation. Each box was vibrated at the chosen 
frequency for 15 minutes. The vibration was resonant and was felt 
through the metallic frame. Following the vibration, - the box was 
sectioned and studied for stress distribution as previously carried 
out in the first phase. 

Tests numbers 3 and 7 ai'e not complete yet; with phases 2 and 3 
of test 3 and phase 3 of test 7 still under Investigation. 

Handling the vibrator was very successful during this period 
of the study and its manual control confirmed the peak that was 
chosen from the pulse analysis for each box. 

Most of the pulse results indicated that a suitable frequency 
within 100 cps can be used to give a peak displacement which reduces 
fatigue of the metal to a minimum level. 
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GAGE H 

1 

2 

3 

4 

LAVER 

1 

0. 

-1265. 

-1265. 

-1265. 

LAYER 

2 

-822. 

-879. 

-879. 

-57, 

LAYER 

3 

-736. 

-791. 

-791. 

-55, 

LAYER 

4 

-650. 

-704. 

-164. 

487, 

LAYER 

5 

-164. 

-1545. 

-1489. 

-1325 

LAYER 

6 

-160. 

-1009. 

-212. 

689 

LAYER 

7 

571. 

-316. 

-206. 

1459 

LAYER 

8 

1970. 

724. 

832. 

list 

LAYER 

? 

688. 

424. 

528. 

1506 


LAVER 10 


U57. 610. 712. 1070 
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LAYER 

J 

6936. 

-3468, -3468. 

6936. 

LAYER 

2 

13965. 

-4674. -4674- 

9348. 

LAYER 

3 

4716. 

-115. -1591. 

6135. 

LAYER 

4 

5231 . 

2276. 

-170. 

5232. 

LAYER 

5 

15160. 

4831. -9823. 

1998B. 

LAYER 

6 

4205. 

-3297.-13716. 

905. 

LAYER 

7 

-4196. 

1700. 

-509. 

7686. 

LAYER 

8 

-2163. 

60. 

.-504, 

9752. 

LAYER 

9 

-1397. 

60. 

-500. 

9473. 

LAYER 

10 

-1514. 

-2041. -2595, 

9898. 

LAYER 

u 

-5008. 

-3686. -4234. 

9083. 

LAYER 

12 

-2230. 

5036. -3231. 

12221. 

LAYER 

13 

2412. 

6386. 

3549. 

8516. 

LAYER 

H 

2691 . 

7231. 

4119. 

9341 , 

LAYER 

15 

412B. 

11548. 10736. 

9908. 

LAYER 

16 

447. 

3557. 

2754. 

5572. 
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LAYER ETCHED 





























CaGi: hi 2 3 4 5 h 7 

LAi'Il^ 1 25737. -19303. '-19303. 28955. 25737. 22520. 25737. 

L,.VEi( 2 13744. -8M7. -10308. 15462- 13744. 12026. 13744. 

LAYER 3 12397. -7805. -9298. 13947. 10961. 9412. 10961. 

LAYER 4 13785. -8672. -8785. 15508. 12174. 12005. 13728. 

LAYER 5 14108,-10412.-10524. 15872. 13995. 10749. 12513. 

LAYER 6 15330. -8991.-10973. 21223. 17090. 14939. 17089. 

LAYER 7 7633. -4588. -4753. 8798. 9023. 5806. 7579. 

LAYER 8 12179. -7677. -7838. 13826. 12241. 10533. 12127. 

-AYER 9 11035. -6547. -6705. 13777. 11095. 8294. 9881. 

LAYER 10 14817.-10298. -8859. 15460. 18064. 12583. 16303. 

LAYER 11 108B1. -6876. -6110. 11511. 11053. 9334. 10832, 

.AYER 12 8611.-5226. -5262. 9836. 9392. 7890. 9790. 

-AYER 13 12713. -7561. -8756. 14524. 12933. 14330. 13940. 

-AYER 14 8412. -5433. -6114. 10283. 10398. 7252. 9125. 
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- y- O.lQOOOOOE-02, • 
y -0.4093776 y 

y 0.1420565 y 
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y 0.5346471 y 
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-y -0.2492783E-01,- 
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Simffliary Of Current In^estli^tlon ‘ 

The research program currently nndervi’ay imrestlgates the 
patterns of residual stresses in welded A3^ steel structures 
of various geometrical designs, formed of hot rolled bars 
of various widths and thicknesses, .High tensile stresses of 
about 25,000 to 30,000 psi were found at the weldi these dropped 
rapidly to compressive stresses in the range of 3,000 to 5,000 
psi which spread for some depth through the bar thickness. 

The pattern of stresses at a distance from the weld was 
reversed, showing high compressive stresses at the outer 
layers. The outer layer stresses, however, were of a lesser 
degree than the stresses found at the weld. 

In the second phase of the study, a load cell was used 
together with an accelerometer to generate an impulse and 
analyze the response, A vibration frequency within the range 
of 100 cps‘ was chosen to generate a peak displacement. Such 
a limited frequency was chosen to reduce the damage to the 
steel by fatigue, due to the vibration. A set of comuuter 
programs was designed for choosing the appropriate frequency. 

This technique was found to be less expensive than using a 
commercial spectrum or Fourier analyzer. Results showed that 
the most effective pulse varies between 50 and 80 cps. 
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At present, the steel structures are vihrated using an M3 
. model CIO mechanical vibrator with attached panel to control 
the frequency to the required value. 

Future Work Plan 

A new study is suisrgested using a variable speed motor with 

( 32 ) 

eccentric load, The motor will be connected to a metallic 
frame fired firmly to the ground, The motor speed can be 
adjusted to generate vibrations at the specified- frequency , 

The study can be performed in the field, or T-rhers the frame 
is to be located. This technique presents a practical way 
to test the method away from laboratory facilities. The 
frequency should be determined using the same present load 
cell and the accelerometer. The pulse can be recorded on the 
screen of a portable oscilloscope and analyzed using the same 
computer programs already developed. 

Fatigue testing of specimens prepared from the sides of the 
welded bores showed a life span of a fatigue limit of 30*000 
psi. This investigation needs to be continued. Similar 
studies are being carried out at present on specimens cut from 
boxes after vibration. Most of the vibrations so far -are 
limited to a maximum of 80 cps-. For the future investigations, 
similar boxes will be vibrated at higher frequencies chosen 

according to other peak values that were registered from the 
impulse analyses. 
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This study '.rill examine the effect of high Tihrations on 
the fatigue life of the steel and establish. a relationship 
between the magnitude of the vibration and the fatigue life 
of the aetal. 

The study will be expanded to include altsminuDi which is used 
on a large scale in the aerospace industry, Alujninirni can 

be etched using caustic soda solution, a tschJiioue that was 

, . (22) 

applied in studying rssiaual stresses in dravjn tubet'. The 

study trill use the same basic eouipsent, including the strain 
indicator, the resistance box and either the existing labo- 
ratory sechanical vibrator or the neu planned variable speed 
motor, lasic and fundamental theories and conp'utsr programs 


presently developed ■will be used 
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This project heretofore has examined the effectiveness of' 
vibration in reducing residual stresses from welding various 
shapes of steel structures. The present activit7, part of 
which is reported here, examined the -oattern of stresses in 
some more steel structures, aluminum structures and the 
effectiveness of vibration in reducing these stresses. 

A second part of the activity studied 'the fatigue life of 
steel exposed to vibration treatment. Further, and as 
reported before, a variable speed motor ha.s been added to 
our testing equipment and was used to vibrate some metallic 
frames at the required frequency. 

A new set of steel frames was formed by weldin<r bars 
h and 8 inch wide by and inch thick, details of shares 
and dimensions are given in Fig. 57. 

The results of stress distributions , mea^sured as before, 
are presented in Figs. 58 -65, , 76 '-82' and 92-96. • ' "The 

second phase covered the load pulse and the acceleration 
response analyses and their results are given in Figs. 66-67. 
83-84 and 97. ■ ■ The three frames were vibrated at 78, 63 

and 112 cps for 15 minutes respectively . 

It was noticed that when the frame was an open structure, 
test2, Figs, 76-82, the peaks of the stresses were lower 
than when the frame was restricted and had a closed, shane, 
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The cross connection between the bars in test 3 helped to 
/ 

reduce the stresses and the results showed a xnaximuBi stress 

of about 5»000 psi only, Pi.‘?s..92-96, 

The alUTuinua (606it 6) studies started by testins? a bar as 

received and showed that t'lc material was exnosed to some 

residual stresses of abo'ab 3'»000 to 4,000 psi which were 

mostly developed due to the straightening process followinec 

the extrusion and treaunent of the bars. The few tests 

carried out so far show that the welding caused peak residual 

3 

stresses of about© 1o fO psi that appeared at a certain 
depth from the outer surface. That may be due to the higher 
conductivity of a,luminuin to the heat of welding as compared 
to mild steel and such a behavior may reduce the risk of 
surface cracks. These conclusions however cannot be final, 

The experiments on welding are still in their first stage and 
further investigations are needed. 

An investigation of how vibration a,ffects the fatisrue life 
of steel was carried out on steel specimens of the shace shovm 
in Pig. 119. The investigation was essential particularly for 
steel exposed to dynamic loading. 

The specimens were machined locally for both the non- 
vibrated and vibre.ted metal boxes. They were tested using a 

classical rotating beam, constant load fatigue testing machine, 
of the type shovm in the same figure. 
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Each specimen was mounted betvreen two rotatini? spindles. 
Stresses were continuously reversed and the steel was passed 
through complete cycles of flexural stress. The load 
consisted of accurately adjusted weights stacked on a hanger 
and the machine operated continuously at a. speed of 10,000 rpm. 

Stress V7as calculated from 
S = 16 wL/tT 

V7here : - 

S=extreme fiber stress, psi 
w=total load on the specimen in lbs, 

L=jnoment arra,i.e,, distance from end support to load 
point (in inches) 

And D=minimuni diameter of specimen (4 inch) 

A maximum fiber stress of 30»000 psi was assumed for the 
structural steel und.er study, for which experiments showed 
a fatigue life span, or cycles to failure of 2x10^ cycles 
and higher. The micros true cure analysis of the grain showed 
no effect of vibrations, Pigs»l20-123. 

The conclusion as of now is that the vibration of steel 
structures vrithin 100 cps and for a limited time period of 
about 15 minutes does not fatigue the steel while it siPhii- 
ficantly reduces the peaks of residual stresses at the weld ' 

zone ; 

Our school of engineering has received a new testing machine 
MTS-Type8lO and I intend to use it in the near future for 
some further studies on the fatigue life of the vibrated 
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metallic structures, both steel and aluminum. 

As reported before, a variable speed motor with a speed 
ran^je of 700 to ^500 has been added to oul* testing 
equipment. It was mounted on a heavy channel-sections base, 
and provided with an eccentric at the end of its spindle. 

The eccentric was connected to the side of the frame to 
be vibrated, Pig. 5.24', The vibration lasted for about 15 
minutes at the required frequency dlcided upon by the impulse 
study and the analysis of the signal registered on the 
oscilloscope. This technique oroved successful and o^jens 
the way for a practical method of vibrating the metallic 
structures at the chosen frequency, in the field as well as 
in the laboratory. 

In conclusion I can say that although the surface stresses 
due to welding mild steel structures are high, yet the 
design of the connections may limit the severity of these 
stresses. The technique of vibrating the structure at 
resonance whithin about 100 cps for a short period of about 
15 minutes reduces the peak values of stresses in both steel 
and aluminum structures. 

As far as fatigue life of the metal, there is no indication 
at the moment of any sound effect of a limited vibration on the 
fatigue of steel, .especially if such vibration is performed 
within 100 cps. 
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Suia-tnary of current investigations 

'The research project, has examined the effectiveness of 
vibration in reducing residual stresses from welding va-rlous 
shapes of steel structures.CL-shaoe connections),. It is 
currently investigating such patterns in welded 606116 aluminum, 
that had been cold extruded. 

Zdigh tensile stresses of about 25,000 to 30,000 psi were 
found at the weld surface of. the steel bars which dropped 
rapidly to compressive stresses in the range of 3,000 to 
5,000 psi that spread for some depth through the bar thickness. 

A cross connection between the bars(T connections) helped to 
reduce the stresses and the results showed maximum surface 
stresses of about 5,000 psi only. It was noticed also that 

when the framis ws.s an open structure, the peaks of the stresses 
were lower than when the frame was restricted and had a 
closed shape. 

Patterns of stresses at a short distance from the vxeld 
were reversed, showing high compressive stresses at the outer 
layers, but to a lesser degree than those foxuid at the weld, 

A load cell and an accelerometer were used to generate an 
impulse in the steel frame and analyze the response picked on 
the screen of an oscilloscope, A vibration frequency within 
the range of 100 cps was chosen to generate a peak disoiacement. 
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Such a limited 'frequency was chosen to reduce the damage; to 
the steel by fatigue , due to vibration. A set of computer 
programs was designed for choosing the appropriate frequency. 

Results showed that the most effective pulse varied 
between 50 and 30 cps. 

Originally the steel structures were vibrated using an 
?I3 model CIO mechanical vibrator x-rith attached panel to 
control the frequency to the required cps* Later a variable 
speed motor with eccentric load was used and connected to 
the metallio frame, Ihe motor speed could be adjusted to 
generate vibrations at the specified frequency, 'This new 
approach was to investigate the possibility of vibrating the 
metallic frames in the field. The pulse was recorded on 
the screen of a portable oscilloscope and analyzed using the 
same computer programs developed previously. 

Patigue testing of specimens prepared from the sides of 
the welded steel boxes were essential particularly for x^elded 
steel structures exposed to dynamic loading. Results showed 
a life span of a fatigue limit of 30,000 psi. The specimens 
for these tests were machined locs,lly for both the non- 
vibrated and the. vibrated metal boxes. They were tested 
using a classical rotating beam . constant load fatigue testing 
machine, where the specimen was mounted between two rotating 
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spindles’ and stresses were continuously reversed. This 
investi.^^ation needs to be continued and a more advanced 
testing; aa.chine MTS-TyueSlO that our department has purchased 
recently will be used. It is also planned to widen the 
study and examine the effect of higher vibrations on the 
fatigue life of the steel and establish a relationship 
between the magnitude of the vibration and the fatigue life 
of the metal. 

The second phase of the project which started few months 
ago deals with the above mentioned aluminum which is being 
used on a large sca.le in the aerospace industry. T-Jhile 
the stress distribution for the steel wa,s studied by etching 
the metal in diluted nitric acid solution, the alumimm is 
being etched using caustic soda solution. Both techniques 
were used before studying the distribution of residual stresses 
in tubes and plates(^^»^^^^5) ^ The study used the same 
equipment used for steel including electric strain gages, 
strain indicator a.nd resistance box and similar computer 
programs modified to suit the aluminum. 

Preliminary studies showed that the as received (cold 
extruded) aluminum bars were exposed to some residual stresses 
of about 3 to ^xlO^ psi which were mostly developed due ’to 
the straightening process following the extrusion and 
treatment of the b3.rs. The few tests carried out so far 
showed peak residual stresses of about 8 to 10x10^ psi that 
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appeared at a certa.in depth from the outer stirfa-ce at the weldi 
•This may be due to the higher conduct ivit 3 r of alumin'um to 
the heat of welding as compared to the mild steel and such 
a behavior may reduce the risj^ of surfa.ce cfacVs, Such 
studies a,re still in their first stage and the conclusions, 
cannot be final. 

Future -'^orh: Plan 

To continue the investigations on aluminum structures. 

To expand the fatigue testing analysis for both the steel 
and aluminum, using the nevr testing machine MTS-Type 810, 

In order to support the findings of the patterns of stress 
distributions, it is necessary to continue the study by 
measuring temperature distribution that takes place -t^hile 
welding and to correlate the patterns of both the stress 
and temperature. Thermoelectric technique using thermo- 
couple ■wires of 0,0005 inch introduced into holes drilled 
in the metallic bars or spot welded to the bar sides. 
Temperatures can be measured at successively decreasing 
distances from the outer surface, together -with temperatures 
at various distances from the weld zone. The state of affair 
may nec.essitate carrying out extrapolation to find the '■ 
temperature at the weld,. In addition a_ radiation pyrometer 
will be used to check on the temperature at the weld. 
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This section of the project studied the effectiveness of 
vibration on releasing residual stresses generated by welding 


aluminum bars to various shapes. The frames were formed of 

1 

if’ 


1 3 1 

bars 3 and 4 inch wide by « and inch thiclc, 'Details 

4 o 2 

are given in Pig. 125 

Results of the stress distribution were worked out as 
before, and are presented in Pi.gures 126-130, 132-141, 143- 
152,. 154-163 and I 65 -I 69 . 


The load pulse and acceleration response analyses are given 
in Figs. 131 ♦ 142, 153 and l64. 

The four frames were vibrated at 82, 69 , 85 and 50 ops 
respectively for 15 minutes. An unsuccessful trial of vibrating 
the frames at 26, 24, 35 and 27 cps, 

was performed for 15 minutes but the embeded stresses were 
not effectively reduced. It was found necessary to vibrate 
at the higher frequencies given above. 

It was observed that the thicker the bars and the more 
restrained the structure (Tests 2&4) the higher the magnitude 
of the trapped stresses. 

Both aluminum and mild steel structures showed comuaratively 
high tensile residual stresses generated by welding. The 
main difference is that peak stresses appeared at the outer 
layers of the mild steel bars, while it appeared at a slight 
depth from the outer layer of the- alTMinum bars. Such a 
behavior is mostly due to the higher conductivity of the aluminum. 
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Previous experimental investigations proved that a mechanical 
vibration of a suitable frequency within 100 cps can be used to. 
reduce the peaks of residual stresses at the weld zone for both 
structural steel and aluminxm frames. The vibration period is 

s ' 

limited to 15 minutes in order to eliminate any fatic^ue in the 
metal bars. 

I 

The earlier experiments with vibration techniques pointed up 

the need to exajSine temperatures at the weld point as they affect 

re''sidual stresses and their distribution. It is believed that 

controlling sharp changes in temperature will modifiy or 

eliminate the harmful effects of residual stresses. 

The basis for the calculations of heat flow originated in the 

work of Robert:!^^^ and Hosentha!^ Various assumptions were 

necessary in the derivations of the mathematical formulae, 

f 1 8 ^ 

including the possible radiation and convection losses' 

The experimental work involved measuring the temperature along 
and across the weld and through 'the denth following the bead 
deposit. Chrome 1-alumel thermocouples '• with refrasil insulation 
were used for that purpose. The Junction of each thermocouple 
was mounted in a hole drilled in the bar to be welded at specific 
locations, as given in the figures with the experimental results. 
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Sometimes the thermocouples were dipped, into the molten metal 

■S’ 

to check the actual melting temperature. Signals from the 
thermocouples were recorded using a pyrometer and a multi-channel 
switch. Equipment has been calibrated by the Omega manufacturers 
and rechecked locally on campus. 

Results are presented as curves of temperature against time. 

The curves have a wide application for different .speeds of moving 
heat source and for differing bar width. Results showed the max, 
temperature alv 7 ays located under the source of heat_ while it was 
moving across the metallic bar. Temperature increased through 
the bar thickness with the successive beads being deposited. 

At the beginning of the weld and during the first pass, a 
steep temperature gradient was always encountered. But as the 
bars got hot the variation of the temperature narrowed. The 
multi bead deposition was carried out on both the A36 steel and 
the 6 o6it 6 aluminum bars using the arc welding technique. 

As the heat source crossed the metallic bar, the temperature 
dropped sharply from the weld zone along the bar, i,e,, in a 
direction at right angles to the weld. This sham drop was due 
to several reasons including the metal conductivity, the sueed 
of weld, the bar proportions and the radiation and convection 
heat losses. 

*A 36 steel indicated a melting temperature of 2300^P 
6061 t 6 aluminum indicated a meltin.g temperature of 1100 °F 
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Figures 172,173- show the variations in temperature in the 
direction of weld and that these variations narrow with the 
successive beads. This supports the idea of limiting the devel- 
opment of stresses if the metal is preheated. It also indicates 
that the larger the number of passes the better the temperature 
distribution and the less the danger of the internal stresses. 

The same line of discussion Is true fox other tests where 
temperatures are measured along the bar , i.e,, at right angles 
to the weld direction, the case that represents the most citical 
problem as far as the variation in temperatures and correspondingly 
the peaks of the generated residual stresses are concerned. 

At present several tests are being ca,rried out on both steel 
and aluminum bars of various thicknesses and widths. This is 
done in order to confirm the temperature distribution character- 
istics and to find out the effect of the bar proportions on the 
degrees of heat losses and the sudden differential temperature. 

For future work, a study will be added to investigate the 
beneficial effects of stress-relief on the metallurgical 
structure ,of both the aluminum and steel and how can this be 
related to the fatigue strength of the metal. An investigation 
will be carried out on the possible reduction or control -of the 
effects of peak temperatures at the weld zone by raising the. 
temperature of the structure members around this zone using some 
type of induction heating thereby reducing the sudden dlffemetial 
of temperature in the surrounding metal. 
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A36 Steel Bar 4 x - 

Gauge arrangement 2 

Holes along the axis o£ the bar 

i\.rc weld m 

Amp. 160, Volt 25, Dia. of rod 1/8 
Speed of weld = -S.S inch/min. 
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During the current experimental investigations, it was observed 
that the thicker the bar and the more restrained the structure, 
the higher the magnitude of the trapped stresses ^ » 3^ ) ^ Peak 

stresses of 25 to 30 x 10^ psi appeared at the outer layers of the 

*7 

mild steel bars, and stresses of 8 to 10 x 10 psi appeared at a slight depth 
from the outer layers of the aluminum bars. Such a behavior is 
attributed to the higher thermal conductivity of the aluminum as compared 
to the steel. 

Experimental analyses proved the effectiveness of the mechanical 

vibration technique on reducing and levelling the residual stresses 

generated by welding in both the steel and the aluminum. The 

vibration period was limited to 15 minutes at a frequency between 

50 and 80 cps, in order to minimize any fatigue of the metal structure. 

Limited mechanical and metallurgical study of the resultant fatigue points 

to further e:^erimental investigations. 

Experimental studies with vibrations pointed out the need to examine 

the temperature at the weld and in the surrounding zone, and how 

the distribution of the temperature affects residual stresses. 

It is believed that by controlling the sharp changes in the temperature 

around the weld it is possible to modify or eliminate the harmful effects 

of the residual stresses. Manual arc welds 3 to 8 inch long are made 

between pairs of steel or aluminum bars - to - inch thick. The bar edges are 

4 2 

straight with 45° chamfer. This procedure represents the normal manual 
arc welding practice. 
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Examination of tenperature distributions both along and across the weld direction, 
is ongoing, using chromel-altmel thermocouples with refrasil insulation. The 
junction of each thermocouple is mounted in a hole drilled in the bar at a specified 
location. Several thermocouples are used simultaneously during each test. Results 
showed that the maximum temperature is located imder the source of heat and moves 
with it across the metallic bar. 

A beneficial finding is that the tenperature residue in the metal increases through 
the bar thickness as the successive beads are deposited. During the initial pass, ^ 
a steep temperature gradient is encountered. As the bar heats the variation of the 
temperature narrows';!. 

Reducing the current in arc welding causes slower melting and shallower penetration. 
A larger number of passes are needed to fill beti\'een the rods, helping to preheat 
the metal and reducing the contrast of temperature distribution. The resulting 
reduction of internal stresses minimizes damage to the metal bars. 
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FUTURE l\iORK PLAN 


At present several tests are being carried out on both steel and aluminum 
bars of various thicknesses and widths. This is being done to cohfiim the 
temperature distribution characteristics and to determine the effect of bar 
proportions on the degree of heat loss and temperature differential. A 
theoretical study of the ten^erature distribution also is tinder way. The 
basis for the calculations of heat flow originated in the work of Roberts 

( 17 ) 

and Rosenthal- • Various assumptions were necessary in the derivation of 
the mathematical formulae, including -the possible radiation and convection 
losses'^®’ 

Mechanical fatigue testing is continuing using -the new MTS-810 "type testing 
machine that has been added to our facilities . This testing mil be 
supplemented with a metallurgical investigation of the grain structure to 
determine -the internal effects of -the vibration treatment. 

A study is needed to analyse the beneficial effects of stress-relief on the 
metallurgical structure of both metals under investigation and how it can be 
related to -the fati^e strength of -the metal. 

It was found that preheating steel to about 400°C produced a reduction in 
residual stresses of about fif-ty percent. The effectiveness of heat trea-fcment 

_( 38 , 39 ) 

following the welding fell off markedly at temperatures lower than 600 C 


Authors findings are based on mechanical tests ra-ther than on actual measurements 

t 

of residual stresses . Future investigation will examine the possible reduction 
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or control of the effects of peak temperatures at the weld zone hy preheating 
the bars around the zone using induction heating. ■ The effect of combining 
such preheating with a limited vibration for about 15 minutes at a frequency of 
50 to 80 cps, following the welding will also be studied. Mechanical and 
metallurgical investigation of the effects of differing vibration frequencies 
on the fatigue life of both steel and aluminum will reveal the optimum range 
for effectively reducing internal stresses. 
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TYPE NMTST.FOR 


00100 


COMPLEX FILL 

00200 


IHTEGER «ORF(9) 

00300 


DIMENSION DH«(15),EQ(15,10) 

00400 


yRITE(5,38) 

00500 

38 

FORHAT(//ENTER H,L,N,G',/) 

00600 


READ(5,*)H,L,N,G 

00700 


URITE(5,599)H 

OOGOO 

599 

FORMAT(^**",F5.0) 

00900 


WRITE(5,39) 

01000 

39 

FORMAT (//ENTER DATA^V) 

01100 


DO 200 1=1, L 

01200 


URITE(5, 201)1 

01400 

201 

F0RMAT(//LAYER',I2) 

01500 


READ(5,*)DHH(I),(E0(I,J),J=1,fn 

01600 

200 

CONTINUE 

01700 


URITE(5,466) 

01800 

466 

FORHAT(///ENTER THE NAME FOR A NEU FILE FOR STORAGES 

01900 


READ(5,112)U0RF 

02000 

112 

F0RHAT(9A1) 

02100 


ENC0DE(9,43,FILL)(U0RF(I),I=t,9) 

02200 

43 

F0RHAT{9A1) 

02300 

02400 

) 

02500 


OPEN(UNIT=20,FILE=FILL) 

«RITE(20,*)H,L,N,G,(LHH(I),I=1,L),((EOa,J),J=1,N),I=1,L 

453 

CL0SE(UNIT=20,FILE=FILL) 

02600 


END 
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TYPE LOOK. FOR 

COMPLEX FILE 
INTEGER UQRK(9) 

DIMENSION E0(15,10),HH(15),BHf1(15),HHn5) 
yRITE(5,l) 

1 F0RHftT(1H3/ENTER THE DATA FILE NA«EO 

READ(5,9?) UORK 
99 F0RNAT(9A1) 

ENC0D£(9,8,FILE)(U0RK(I),I=1,9) 

8 FORHAT{?A1) 

OPEN(UNIT=20,FILE=FILE) 

89 READ(20,*,END=90) H,L,H,G, (BHM(I) , 1=1 ,L) , ( (E0(I, J) , J=1 , 

1 N),I=1,L) 

90 HHH=0.0 

DO 50 1=1, L 
HHa)=HHM+DHM(I) 

HHN=HH{I) 

HH(I)=H-HH(I) 

50 CONTINUE 

888 CL0SE(UNIT=20,FILE=FILE) 

URITE(5,21)H,N,L,G 

21 F0RHAT(//,5X,' BAR THICKNESS%F6.0,5X,"8 

1 OF GAGES^I3,5X/^1 OF LAYERS" , 13, 5X/GAGE SE 
1TTIN6",F8.2) 

URITE(5,22) 

22 F0RHAT(//,5X, "AVERAGE THICKNESS", 5X, "AVERAGE 

1 REMOVED", 5X, "GAGE NUMBER 1 2 3 

1 4 5 6 7") 

DO 45 11=1, L 

HRITE{5,26)HH(I1),DHM(I1),(E0(11,J),J=1,N). 

26 F0RMAT(/,10X,F5.0,18X,F4.0,23X,10F8.0) 

45 CONTINUE 

END 



-TYPE 

0SI00 

03200 



03600 
03700 
00800 
00900 
01000 
0!100 
01200 
01300 
01400 
01500 
01600 
01700 
01800 
01900 
02000 
02050 
02060 
02100 
02200 
02300 
02400 
02500 
02600 
02700 
' 02800 
^ 02900 
03000 
03100 
03120 
' 03200 
03300 
03400 
03500 
03630 
03700 
03900 
04000 
04100 
04200 
04300 
04503 
04600 
04700 
04800 
04900 
05000 
05100 
05200 



15500 

"05600 

05700 

05800 

05900 

06000 


M\iE 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


H IS THE ORIGINAL THICKNESS (IN.) 

HH IS THE DEPTH OF THE I LAYER REHOVEB 
DHH(I) IS THE THICKNESS OF THE I LAYER REHOyEEi 
EO(I,K) IS THE CHANGE IN STRAIN IN THE OUTSIDE FIBER 
I IS THE LAYER 
K IS THE LOCAIION 

STRSS(I,K) IS THE TOTAL STRESS IN THE LAYER REHOUED 
I IS THE LAYER 
K IS THE STRAIN LOCATION 
SSUN(I,K) IS THE SUN OF THE SUPERIHPOSEB STRESS 
I IS THE LAYER 
K IS THE STRAIN LOCATION 
SRES(I,K)- IS THE RESIDUAL STRESS IN A LAYER 
I IS THE LAYER 
K IS THE LOCATION 
L IS THE NUMBER OF LAYERS 
N IS THE NUMBER OF STRAIN LOCATIONS 
HHfl IS AN INTERNAL COUNTING ','ARIABLE 
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COMPLEX FILL 
INTEGER UORF(?) 

COMPLEX FILE 

DIMENSION HM(19),IiHMn9),HN(19),DHN(19),EOn9.10),STR3S( 
1 19,10J,SSUM(19,10),SRES(19,10),HH(19) 

URITE(5,1 ) 


1 

FQRHAT(1H3,' ENTER THE DATA 
CALL OPEN(20,IER,FILE) 
IFdER .EQ. 0)60 TO 87 
60 TO 89 

FILE NAME") 

CO 

URITE(5,8B) 


'88 

FORMAT!/, ' ERROR") 


89 

REAIII20, (•,ENH=495) H,L,N,G, 

!DHH!I),I=1,L),!!E0II 


1,J),J=1,N),I=1,L) 

— 

495 

DO 15 J=1,L 
DO 10 K=1 ,N 


10 

SSUH!J,1!)=3.0 


■ 15 

■ CONTINUE 
HHM=0.0 
DO 50 1=1, L 



HM!I)-=HHM+IlHH!n ’ 



HN!I)=HH!I) 

DHN!I)-=HHri!I) 

HHH=HM!I) 



HH!I)-=H-HM!I) 


50 

CONTINUE 



CLOSE (UNIT=23,FILE=FILE) 
no 500 1=1, L 
DO 700 K=1,N 

STRSS(I,K)=29.0-)EO(I,K) l (~1 .0}/(DHM(I)/(H-HN(I) )-3.0-» 
IDHH(I-) KH-HHm+riHH(I) )/(H-HM(I))-i 12) 

DO 600 J=I+1,L 

600 SSUH{J,K)=29.0+EO(I,K) M (H-HM(D)-i 12+3.0 MH-Hm I ) 

1+riHH(I))MH-2.0iHN(J)+DHN(J) + HHa)))/((H-HM(I))-Ji2~3.0i(H- 
1HHm + HHHa))i (H-HM(I)))+S3UM(J,K} 

700 SRES(I,K)=STRSS(I,K)-3SUli(I,l() 

500 CONTINUE 

URITE(5,21 )H,N,L,6 

21 FORMAT!//, 5X," BAR THICKNESS" ,F6 .0 ,5X, '8 

1 OF GAGES', 13, 5X, "It OF LAYERS ', 13 . 5‘‘ 'GAOE 3E 
1TTING',F8.2) 



T T 

26 

45 

66 

34 

47 

898 

466 

112 

43 


1 

99 

999 

8 

9 

10 
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fORf1AT(//,5X, ' AyilRAGE 1HICKNE33 ' ,5X , AVERAGE RE 
1 MOVED ',5X, 'GAGE HUMBER 1 2 3 

14- 5 6 7 ) 

DO 45 11=1. L 

URI f E ( 5 . 26 ) HH ( 1 1 ) , HH,4 ( 1 1 ) , ( EO < 1 1 , J ) , J~ 1 , ^ 
FORmAK/, ]0X,F5.0, 18X,F4.0,23X,13F8.0) 

CONTINUE 

URITE(5,66) 

FORMAT (////, IH:/ gage 11 1 2 3 

1 5 6 7 8^) 

DO 47 I1=1,L 

URITE(5,34m ,(SF:ES(I1 ,K),K=1 ,N) 

FORMATC/, ' LAYER-, 3X, 12, 3X,10F7.!3) 

CONTINUE 
TYPE 893 
FORMA rn M2) 

UftITE(5,466) 

FORMAT (//,•' ENTER THE NAME FOR A NEU FILE FOR 
1STORAGE--) 

REAIK5, 1 12)U0RF 
FORMAr(9AI) 

EHC(IUE(9,43,FILLHU0RF(I),I=1,9) 

FORMAT(9A1) 

OPEN(UNIT=20,FILE=FILL) 

URITE(20, i)H,L,N,G,(fiHN(I) ,I = 1,U,{ (E0(I,J),J=1 , 
1N),I = 1 ,L), ([1H(I),I = 1',L) , nSRES(I,J),J = 1 ,N),I = 1 ,L) 
CLOSE(-UNIT=20,FILE=FILL) 

END 

SUBROUTINE DPENIUNIT, IER,FILE) 

INTEGER UNIT 
INTEGER ERROR 
INTEGER IJ0RK(9) 

COMPLEX FILE 
INTEGER COL 
DO 1 C0L=1,9 
UORK(COL)-=' 

CONTINUE 

. COL=1 

READ(5,9?)U0RK 
F0RUMAT(?A1 ) 

EliC0DE{9,8,FILE)(UQfii;m ,I--1, 2) 

FORMAT! ' -,2A5) 

FORMAT CyAT) 

0P£N(UNIT=UN!T,FILE=FILE,ERR--9) 

IER = I 

RETURN 

END 
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TYPE HAKA.FOR 

H IS THE ORIGINAL THICKNESSdN. ) 

HH IS THE DEPTH OF THE I LAYER REHOMED 
DHM(I) IS THE THICKNESS OF THE I LAYER REMOVED 
E0(I,K) IS THE CHANGE IN STRAIN IN THE OUTSIBE FIBER 
I IS THE LAYER 
K IS THE LOCATION 

STRESS(I,K) IS THE TOTAL STRESS IN THE LAYER REMOVED 


I IS THE LAYER 
K IS THE STRAIN LOCATION 
SSUH(I,K) IS THE SUM OF THE SUPERIMPOSED STRESS 
I IS THE LAYER 
K IS THE STRAIN LOCATION 
SRES(I,K) IS THE RESIDUAL STRESS IH A LAYER 
I IS THE LAYER 

K IS THE LOCATION 
L IS THE NUMBER OF LAYERS 
N IS THE NUMBER OF STRAIN LOCATIONS 
HHH IS AN INTERNAL COUNTING VARIABLE 

COMPLEX FILL 
INTEGER «0RF(9) 

COMPLEX FILE 

DIMENSION HM(19),DHM(19),HN(19),DHN(19),E0(19,1O),3TRSS( 
119,10),SSUM(19,10),SRES(19,10),HH(19) 
yRITE(5,1) , 

i FORHAKIHS," ENTER THE DATA FILE NAME') 

CALL 0PEH(20,IER,FILE> 

IFdER ,EQ. 0)G0 TO 87 
GO TO 89 

87 URITE(5,88) 

88 FORHATt/," ERRORS 

89 READ(20,*,ENDM95) H,L,N,G , (DHM( I) , 1=1 ,L) , ( (E0( I 

1,J),J=1,N),I=1,L) 

495 DO 15 J=1,L 

DO 10 K=1,H 

10 SSUM{J,K)=0.0 

15 CONTINUE 

HHM=0.0 
DO 50 1=1 ,L 
HM(I)=HHM+DHH(I) 

HN(I)=HM(I) 

DHN(I)=DHM(I) 

HHH=HH(I) 

HH(I)=H-HH(I) 

50 CONTINUE 

CL0SE(UNIT=20,FILE=FILE) 

DO 500 1=1, L 
DO 700 K=1,N 




STRSS(I,K)=10.0t=EG(I,K)*(-1.0)/(DHH(I)/(H-HH(I))-3.0:*= 

DO 600 J=I+1,L 

' 600 SSUH(J,K}=tO.O*EOa,K)i((H-HH(I))**2+3.0t(H-HH(I) 

1+DHHa))*(H-2.0*HN(J)+DHN(J)+Hiia)))/((H-HH(I))+*2-3.0*(H- 

1HH(I)+DHH{I))*(H-HH(I)))+SSUH(J,K) 

700 SRES(I,K)=STRSS(I,K)-SSUH(I,K) 

500 CONTINUE 

gRITE(5,21)H,N,L,B 
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21 F0RHAT(//,5X,"BftR THICKNESS', F6.0,5X,"# 

"■ '1 OF GAGES", 13, 5X,"S OF LAYERS"',I3",5X",^GSGE SE 

1TTINB",F8.2) 

«RITE(5,22) 

22 F0RHAT(//,5X, "AVERAGE TH I CKKESS",5X," AVERAGE RE 

1M0VED",5X,"6AGE KUHBER 1 2 3 

1 4 5 6 7") 

DO 45 11=1, L 

HRITE(5,26)HHai),DHN(I1),(E0(I1,J)-,J=1,N) 

26 FDRHAT(/,10X,F5.0,18X,F4.0,23X,10F8.0) 

45 CONTINUE 

URITE(5,66) 

66 F0RHAT(////,1H2," GAGE It 1 2 3 

1 4 5 6 7") 

BO 47 11=1, L 

«RITE(5,34)I1 ,<SRES(I1,K),K=1,N) 

34 FORHATI/, "LAYER", 3X, 12, 3X,10F7.0) 

47 CONTINUE 

TYPE 893 

898 F0RHAT(1H2) 

WRITE(5,466) 

466 FORHATI//, 'ENTER THE NAME FOR A NEU FILE FOR 

1 STORAGE") 

READ(5,112)U0RF 
112 F0RHAT(9A1) 

ENCQDE(9,43,FILL)(U0RF(I),I=1,9) 

43 F0RHAT<9A1) 

0PEN(UNIT=20,FIIE=FILL) 

«RITE(20,=t')H,L,N,G,1DHN(I),I=l,U,((ED(I,J),J=l, 
1N),I=1,L),(HH(I),I=1,L),((SRES(I,J),J=1,N),I=1 ,L) 
CL0SE(UNIT=2d,FILE=FILL) 

END 

SUBROUTINE. OPEN(UNIT,IER,FILE) 

INTEGER UNIT 
INTEGER ERROR - 
INTEGER U0RK(9) 

COHPLEX FILE 
INTEGER COL 
DO 1 COL=1,9 
«QRK(COL)=" " 

1 CONTINUE 

C0L=1 

READ(5,99)«0RK 
99 F0RHAT(9A1) 

ENC0DE(9,8,FILE)(U0RK(I),I=1,9) 

999 FQRHATI" ",2A5) 

8 FQRHAT(9A1) 

. . 0PEN<UNIT=UNIT,FILE=FILE,ERR=9) 

9 IER=1 

10 RETURN 
END 


AOAftOJ‘‘7 
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WAl Main program for frequency response analysis 

MA2 Subroutine for transfer function (complex valued) 

and coherence function (Real). It makes complex 
variables out of real raw and smooth versions. 

VA3 Subroutine for power spectral density functions 

(Raw version and smooth version by Parzen Lag). ■ 

WA4 Subroutine for cross spectral density function 

(Raw version and smooth version by Hanning). 

HA5 Subroutine to check data and average X & Y (made 

equal to zero) 

■!'rA6 Subroutine used to make the input and output files 

(Opens, Reads and Closes the files containing data) 

INTER Intepolates X & Y data 

FRFOUR Finds Fourier Transfer. .SUBMIT FRFOUR Prints 
TRANS (Includes POUR). 

FRUA Finds statistical transfer values (OUTPUT) . 

SUBMIT PRVTA Executes UA Series and the output 
includes* OUTPUT, POR20.DAT, POH30.DAT, 

FOR31.DAT and POR32.DAT (No TRANS. )(M value is needed) 

FRR Finds both Fourier' Transfer Values (TRANS) -and 

the Statistical Transfer Values (OUTPUT) (M value is needed) 


Possible Output Files 
FOR20.DAT Results of averages 

F0R30.DAT X power spectrum results 

F0R31.DAT- Y power spectrum results 

F0R32.DAT Cross spectrum results (Comulex) 
FOR33 CK(I),QK(I) - 

FOR36 . AR(I),BR(I) 

P0R3 ? CKCX ( I ) , QKCX ( I ) Raw 

FOR38 CKCX ( I ) , QKCX ( I ) Smooth 

POR55 CPK(I),QPK(I) 


*Pile Name For Storage Of Results Must Be 9 Characters 
With .DAT Ext. 



TYPE 

00100 

00200 

00300 

00400 

00500 

00600 

00700 

00800 

00900 

01000 

01100 

01200 

01300 

01400 

01500 

01600 

01700 

01S00 

01900 

02000 

02100 

02200 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 


OF 


QUAl-ttY 


UA1 .FOR- 

PROGRAM VIBilAT 

COHHON GYX(100),GXYS(100) 

COHPLEX FILEX,HT,HS,GXY,GXYS 

DIHENSION X(999),Y(999),GKX(100),GKXS(100),GKY(100), 

1 GKYSnoO),HTnOO),HSnOO) jC(100),CS(100) 

1 ,D(100),RXn0O) 

1 ,A(100),BB(100),CC(100),CD(100),CSS(100),DC(100) 

0 ENTER THE LAG NUMBER, M 

TYPE 5 

5 FORHATC' H? ') 

READ(5,*)M 

C OPENS READS, AND CLOSES THE FILE CONTAINING THE 

C DATA FOR THE INPUT FUNCTION, X 

CALL NAHE(FILEX) 

OPEN{UNIT=20,FILE=FILEX) 

REAB(20,*)HX,NX,TR,(X(I>,I=1,NX) 

-40 CL0SE(UNIT=20,FILE=FILEX) 

C CHECKS TO MAKE SURE M, AND N ARE PQSITIOE AND EVEN 

C CALCULATES CUTOFF FREB BANDUIDTH, AND STANDARD ERROR 

CALL CHECK(NX,HX,H,F,B,E) 

»RITE(5,50)NX,HX,M,F,B,E 

URITE<46,50)NX,HX,M,F,B,E 

50 F0RiiAT(5X,''NX='',I3,3X,''H='',F15.5,3X,''M=", 

1 I3,3X,TC=^F15.5,3X,•'BE=^F15.5,3X,''E='',F15.5) 

C READ DATA FOR OUTPUT FUNCTION, Y 

CALL NAME(FILEX) 

OPEN(UNIT=20,FILE=FILEX) 

READ(20,*,END=75)HY,NY,TR,<Y(I),I=1,NY) 

75 CL0SE(UNIT=2O,FILE=FILEX) 

C CHECK Y DATA 

CALL CHECK(NY,HY,H,F,B,E) 

WRITE(5,51)HY 

51 FORMAT!" NY" ,13) 

C HAKE SURE THERE ARE SAME NO. OF X AND Y VALUES 

IF1NX .NE. NY) STOP "NX NOT EQUAL TO NY" 

C AVERAGE X AND Y HADE EQUAL TO ZERO 

C RESULTS OF AVERAGING IN FILE F0R2O.DAT 

CALL FORH(XjNX) 

CALL FORH(Y,NY) 

WRITE(20,291 

29 FORMAT (3X,"X(I)",5X,"Y(I)",5X, "I") 

«RITE(20,31}(X{I),Y(I),I=1,NX) 

31 F0RHAT(5X,I3,10X,2E10.5,F15.5) 

C POUER SPECTRA FOUND FOR X AND Y 

C RESULTS FOR X FUNCT. FOR30.DAT 
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04800 

04900 

05000 

05100 

05200 

C 

05300 

C 

05400 

05500 

05600 

C 

05700 

C 

O5S0O 

05900 

06000 

C 

06100 

C 

06200 

06300 

06400 

06500 

06600 

C 

06700 

06800 

06900 

101 

07000 

07100 

105 


RESULTS FOR' Y FUNCT. ‘ F0R31 7MT 
CALL POUER(X,HX,NX,h,DKX,GKXS,S,RX) 
gRITE( 30 , 4 =)(GKX(I),GKXSa),I= 1 ,H) 

CALL POHER(Y,HY,NY,H,GKY,GKYS,L,RX) 
URITE(31,*)(GKY<I),GKYS(I),I=:1 ,H) 

COHPLEX VALUED CROSS-SPECTRAL DENSITYFOUND 
RESULTS IN FILE F0R32.DAT 
CALL CR0S3(X,Y,HX,N,NX,A,BB,CC,CS,CS3,BC) 
URITE(32,*)(GXYn),GXY3<I),I=1,H) 

TRANSFER FUNCTION FOUND, COMPLEX VALUED 

COHERENCE FUNCTION FOUND. REAL 

CALL TRANS (GKX,GKY,GKXS,QKYS,N,HT,HS,C,C3, 

A,BB,CC,CD) 

FILE IS NAHED FOR STORAGE OF THE RESULTS 
FILENAME MUST BE 9 CHARACTERS UITH .DAT EXT 
TYPE 105 
CALL NAHE(FILEX) 

0PEK(UNIT=20,FILE=FILEX) 

URITE(20,101) 

F0RliAT(5X"HT'',5X,'’HS(I)',5X,^Ca),',5X,'CS<I)'-,5X,aV//) 

URITE(20,*)(HT(I),HSa),C(I),CS(I),I=1,ii) 

CL0SE(UNIT=20,FILE=FILEX) 

FQRMAT("NAHE OF FILE FOR TRANSFER DATA?'') 

END 



TYPE 

OOlOO 

00200 

00300 

00-^00 

00500 

00600 

00700 

00800 

00900 

01000 

01100 

01200 

01300 

01400 

01500 

01600 

01700 

01800 

01900 

02000 

02100 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 

04800 

04900 

05000 

05100 

05200 

05300 

05400 

05500 

05600 

05700 
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UA2.F0R 

SUBROUTINE TRANS(GKX,8KY,6KXS,GKYS,H,H,HS,C,CS, 

1 6KXC,CX,CXS,GKXSC) 

C THIS PROGRAM FINDS THE TRANSFER AND COHERENCE FUNCTION 

C COEFFICIENTS BETUEEN TUG SIGNALS. 

EXTERNAL ABS 

COMHON GXY(100>,GXYS(100) 

COMPLEX GXY,GXYS,GKXC,GKYC,GKXSC,BKYSC,H,HS,CX,CXS 
REAL C,CS 

DIMENSION GKXC(H),GKXS(ti),GKX(H), 

1 GKY(M),GKYS(H),H(H),HS(H),CX(H),CXS(H),C(H),CS(M) 

1 ,GKXSC(H) 


C FIND THE TRANSFER FUNCTION 

N=H+1 

DO 15 K=1,H 

IFJGKX(K) .EQ.O)GKX(K)=.0001 
IF(GKXS(K).EQ.O)GKXS(K)=.0OO1 
IF(GKY(K).EQ.0)GKY(K)=.OOO1 
IF(GKYS(K).EQ.O)GKYS(K)=.0001 
C MAKE COMPLEX VARIABLES OUT OF REAL 

GKXC(K)=GKX(K) 

GKXSC(K)=GKXS(K) 

C THE "RAU VERSION" 

H(K)=GXY(K)/6KXC(I() 

C THE "SMOOTH VERSION" 

HS(K)=6XYS<K)/6KXSC(K) 


C THE COHERENCE FUNCTION 

C THE RAU VERSION 

CX(K)=GXY(K)-s-GXY(K) 

C(K)=C(K)/GKX(K)/GKY{K) 

C(K)=ABS(C(K)) 

C THE SMOOTH VERSION 

CXS(K)=6XYS(K)*6XYS(K) 

CS(K)=CXS(K) • 

CS(K)=CS(K)/6KXS(K)/GKYS(K) 

CS(K)=ABS(CS(K)) 

15 CONTINUE 

RETURN 

END 

CHNittlltMIINIItINNItMMPjHHINililHMIMNNNiHHHMHHIHIHHHI 



TYPE 
13400 
13500 
13600 
13700 
13800 
14000 
14100 
14200 
14300 
14400 
14500 
14600 
14700 
14800 
14900 
'TS'IOO 
15200 
15300 
15400 
15500 
15600 
15700 
15800 
15900 
16000 
16100 
16200 
16300 
16400 
16500 
16600 
16700 
16800 
16900 
17000 
17100 
17200 
17300 
17400 
17500 
17600 
17700 
17800 
17900 
18000 
18100 
18200 
18300 
18400 
18500 
18600 
18700 
18800 
18900 
19000 
19100 
19200 
19300 
19400 
19500 
- 19600 
19700 
19800 

19900 


UA3.F0R 

SUBROUTINE POyER(X,H,N,H,GK,GKS,IlR,RXR) 


DIHENSION X(H),GK<N),eKS(H),RXR(i1),DR(Jn 
EXTERNAL COS 


IlO 20 IR=1,H 
RXR<IR)=0. 

NX=N-IR 
DO 10 H1=1,HX 

RXR(1R)=RXR(IR)+X(N1 )=t=X(N1 + IR) 
10 CONTINUE 

RXRUR)=RXR(IR)/flX 
20 CONTINUE 


puss 


R0=0. 

DO 30 N1=1,N 
RO=RO+X(N1 
30 CONTINUE 
RO=RO/N 


MM=H-1 

QM=ii 

MK=H+1 


C RAW VERSION 

DO 40 K=1,HK 
Q=K-1 

Q=Q/QH={-3.I416 

GK(K)=0. 

DO 50 IR =1,H« 

R=IR 

GK(K)=8K(K}+RXR(IR)=i=CaS(R+a> 

50 CONTINUE 

BK(K)=2.*Hi=(R0+2.+BK{K) + (-1.**(K-n)=i^RXR{H)) 
40 CONTINUE 

HH=ii/2 

Q=H 

C SHOOTH BY PARZEN LAG FOR SHOOTH VERSION 

DO 60 IR=1,HH 
R=IR 

DRCIR)=1 .-6.*{ (R/Q)4:*2-(R/a)i=-i=3) 

60 CONTINUE 

BR(l>=l. 

DO 70 IR=HH+1,H 
R=IR 

DR(IR)=2.>t:((1.-R/Q}=i::l:3) 

■70 CONTINUE ^ 


DO 80 K=1,NK 
Q=K-1 

Q=a/aH*3.1416 

GKS(K)=0. 

DO 90 IR=1,MH 

GKS(K)=GKS(K)+RXRaR)^Dft(IR)tCOS(RisQ) 
90 CONTINUE 

GKS(K)=2.*H=KR0+2.+GKS(K) ) 

80 CONTINUE 


RETURN 

FMIi 



TYPE 

00100 

00200 

00300 

00^00 

00500 

00600 

00700 

00800 

00900 

01000 

01100 

01200 

01300 

01400 

01500 

01600 

01700 

01800 

01900 

02000 

02100 

02200 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 

04800 
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UM.FOR 


subroutine cross (x,y,h,rt,Ti>!,br,cir,ry>:,r>{y,*:ik,cU) 

« 

diHension x(dk) ,y (nx) ,ck(«) ,qk,-(tt> , 

1 RXY<H),RYX(H),AR(H),BR(H},QKCX(10O),CKCX(1CO) 
coHMon gxy(100),5xys(100) 

coftplex gxy,gxy5,gkcx,qkcx,ckcx 

do 10 ir=0,ft 

i=ir+1 

nr=nx“ir 

rxy(i)=0. 

ryx(i)=0. 

do 20 n=1,m’ 

rxy ( i ) -rxy ( i ) ) it*y (m+ir ) 

ry::(i)=ryx{i)+y(n)*>i(n+ir) 

20 continue 

rxy(i)=r:iy(i)/nr 
ryx<i)=ryx(i1/nr 
ar ( i ) =0 . 5i= { rxy ( i ) +r yx ( i ) ) 
br(i)=0.5=t'(rxy{i)-ryx(i) ) 

10 continue 

do 15 k=0,H 

q=k 

qH=« 

q=q/qft*3.1416+.001 

kk=k+1 

ck(kk)=0.0 

qk(kk)=0.0 

do 25 ir=2,H 
r=ir-1 

ek(kk)=ck1kk)+2.=Kar(ir)*cos(q*r) ) 
qk(kk)=qk(kk)+br(ir)i=5in(q+r) 

25 -continue 

CK(KK)=2.*H*(AR(l}+CK(KK) + <-1.-t=fK)-tAR(H+1)) 
qk(kk)=4.*h*qk(kk) 

15 continue 

Hft=rt+1 



original PASe e®’ 

op POOR 


04900 



05000 

05t00 


do 30 i=l,HH 

05200 


qkc>{(i) = (0,-l ) 

OS300 


ckcx(i)=ck(i) 

05400 


qkex(i)=qkc>;(i)*qk(i) 

05500 


gxy(i)=ckcx(i)+qkcx(i) 

05600 

05700 

30 

continue 

05800 

c 

SMOothing by banning oethod 

05900 


ck(l)=(ck(l)+ck(2))/2. 

06000 


qkd )=(qk(l)+qk(2))/2. 

06100 


do 37 j=2,ft 

06200 


ck(j+1)=(ck(j)+ck(j+1))/2. 

06300 


QK(J+l)=(QK(J)+QK(J+1))/2. 

06400 

37 

continue 

06500 


do 40 i-2,rt 

06600 


ck(i)~(ck(i-l)t2.+ck(i)+ck(i+1))/4. 

06700 


qka>=(qk(i-1)+2.=i'qk(i)+qk(i+1))/4. 

06800 

06900 

40 

continue 

07000 

c 

SHOoth version 

07100 

07200 


do 50 i=l ,HH 

07300 


ckcx(i)=ck(i) 

07400 


qkcx(i)=(0,*1 ) 

07500 


qkcx(i)=qkcx(i)i=qk(i) 

07600 


GXYSa)=CKCX{I)+BKCX(I) 

07700 



07800 

07900 

50 

■ continue 

08000 


return 

08100 


end 
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TYPE WAS. FOR 


00100 


SUBROUTINE F0RM(X,N) 

00200 


BIHENSION X(H) 

00300 


SUH=0. 

00400 


SN=N 

00500 


DO 10 1=1, N 

00600 


SUH=SUH+X(I) 

00700 

10 

CONTINUE 

00800 


SUH=SUi1/SN 

00900 


DO 20 1=1, H 

01000 


X(I)=X(I)“SUH 

01100 

20 

CONTINUE 

01200 


RETURN • 

01300 


END 

01400 


01500 

01600 


SUBROUTINE CHECKIN, H,H,FC, BE, E) 

01700 


EXTERNAL ABS,SQRT 

01800 


HH=ABS<H/2) 

01900 


HN=2-(=HH 

02000 


IF(HN .NE. H)PAUSE "M NOT POSITIOE, EVEN'' 

02100 


NQ=ABS(N/2) 

02200 


NC=2*NQ 

02300 


IF(NC .NE. N)PAUSE "N NOT POSITIVE, EVEN'' 

02400 


FC=.4/H 

02500 


QH=H 

02600 


BE=1./(QH=«=H) 

02700 


QN=N 

02800 


E=SQRT(QH/QN) 

02900 


RETURN 

03000 


END 

03100 

C 
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TYPE 

00100 

UA6.F0R 

SUBROUTINE NAHE(FILEX) 

00200 


COHPLEX FILEX 

00300 


INTEGER U0RK{9) 

00400 


TYPE 30 

00500 

30 

FORMAT ("FILE NAilEsDATA;") 

00600 


READ(5,99HW0RK(I),I=1,9) 

00700 

99 

F0RHAT(9A1) 

00800 


ENC0DE(9,8,FILEX)(«0RI<(I),I=1,9) 

00900 

8 

F0RHAT(9A1) 

01000 


RETURN 

01100 


END 

01200 




00100 

00200 

00250 

00300 

00400 

00500 

00600 

00700 

00800 

00900 

00950 

01000 

01100 

01150 

01UO 

01200 

01300 

01400 
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PROGRAM FOUR 

BIHEMSION X(999),Y(999),H(999UA(99?) 

CALL ERXSETdOjO) 

READ(23,=s,END=33) HP,MX,T, (X(I) ,1=1 ,NX) 

33 READ(24,=i=,ENB=44) HP, NY, T, (Yd) ,1=1 , NY) 

44 CALL F0RM(X,NX) 

CALL F0RH(Y,NX) 

CALL F0URTR(X,NX+4,-1 ,A) 

CALL F0URTR(Y,RX+4,-l ,A) 

BO 50 1=1, NX 
IF(Xa).EQ.O.)X(I) = .O01 
Ha)=Ya)/X(i) 

50 Aa)=I/T 

HRnE<41,11) 

11 F0RMAT(4X,'Aa)'‘,10X,^X(I)',1OX,"Y(I)M0X,'H{I)^///) 

URITE(41 ,*)(Aa),X(I),Y(I),H(I),I=1,HX) 

10 FCRHAT(10X,F10.3,5X,3E15.12) 

ENB 
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00100 
00200 
00300 
00400 
00500 
00600 
00700 
00800 
00900 
OlOOO 
01100 
01200 
01300 
"O'l 350 
01400 
01500 
OUOO 
01700 


■PROGRftH INTERS 
DIHENSION Y(999),X(300) 

• READ(23, + )A,J,C,(Xa),I=1,J) 

TYPE 5 

5 FORHATI' INPUT NUMBER OF INTERPOLATED VALUES") 

READ<5,*)N 
L=N*J 
A=A/N 

DO 10 1=1, N 
K=0 

DO 15 JJ=1,L,N 
K=K+1 

Y(JJ+I)=X(K)+I=i=(X(K+1)-X(K))/N 
: ■■■ YUJ)=X(K)' ' " " “ 

15 CONTINUE 

10 CONTINUE 

HRITE(24,*)A,L,C,(Y(I),I=1 ,L) 

END 
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.TYPE FRFOUR.FOR 

COPY F0R23.BAT=X.DftT 

COPY FQR24.BAT=Y.BA7 

EX F0UR.REL,yA5.REL,SYSsSPCLIB.REL 

DELETE TRARS 

RENAME TRARS=F0R41.DAT 

PRINT TRANS 

DEL F0R23.DAT, F0R24.DAT 

DIR/L 

DAYTIME 



OF POOR QUAUTf 


.TYPE FRUA.FOR 

COPY FQR23.BAT=X.BAT 

COPY FDR24.BAT=Y.DAT 

COPY F0R60.HAT=H.»AT 

EX UA11,UA2,UA3,UA4,UA5 

DEL FQR60 . DAT , F0R23 . DAT , F0R24 .BAT 

PRINT F0R20. DAT, F0R30. DAT, F0R31 .DAT, F0R32.DAT 

DEL OUTPUT 

RENAHE- 0UTPUT=F0R46.BAT 
PRINT OUTPUT 

DEL FQR20.DAT,F0R30.DAT,F0R31 .BAT, F0R32.DAT 

DIR/L 

DAYTIHE 
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.TYPE FRR.FOR 

COPY FDR23.DflT=X.BAT 

COPY F0R24.DAT=Y.BAT 

■E)TF0UR.REL,«A5:REL',SYS;'SPCLIB'.REL 

DEL TRANS 

RENAME TRANS=FOR41.5AT 

PRINT TRANS 

COPY FDR60.DAT=H.DAT 

EX UA1t,UA2,UA3,UA4,UA5 

DEL F0R60.DAT, F0R23.BAT, FQR24.DAT 

PRINT F0R20. SAT, F0R30.BAT, F0R31.BAT, FQR32.DAT 

DEL OUTPUT 

RENAME 0UTPUT=F0R46.DAT 
PRINT OUTPUT 

DEL F0R20.DAT, F0R30.DAT, F0R31.DAT, FCR32.BAT 
DIR/L / 

DAYTIME 
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TYPE 

00100 

UUA4.F0R 


00200 


SUBROUTINE CROSS(X,Y,H,M,NX,BR,AR,RYX,RXY,QK 

00250 

00300 

1 

,CPK,QPK) 

00400 


DIHENSIDN X<HX),Y{NX),CK(H),QK.(H), 

00500 

1 

RXY(«),RYX(H),AR(H),ER(H),BKCX(100),CKCX(100) 

00550 

1 

,CPK(100),BPKn00) 

00600 


C0H«0N GXY(100),8XYSU00) 

00700 

00800 


COMPLEX 6XY,6XYS,GKCX,0KCX,CKCX 

00900 


DO 10 18=0, H 

01000 


I=IR+1 

01100 


N8=«X-IR 

01200 


RXY(I)=0. 

01300 


RYX(I)=0. 

01400 



01500 


DO 20 N=1,NR 

01600 


RXY(I)=RXYa)+XCM)*Y(N+IR) 

01700 


RYX(I)=RYXa)+Y(H)*X{N+IR} 

01800 

01900 

20 

CONTINUE 

02000 


RXY(I)=RXY(I)/NR 

02100 


RYX(I)=RYX{I>/NR 

02200 


AR(I)=0.5=f=(RXY(I)+RYX(I)) 

02300 

02400 


BRa)=0,5*(RXYa)-RYX(D) 

02500 

10 

CONTINUE 

02550 

02600 


URIT£(36,=k}(AR(I),BR{I),I=1,M+1> 

02700 

02800 


DO 15 K=0,H 

02900 


Q=K 

03000 


QH=H 

03100 


Q=Q/QH*3.1416+.O01 

03200 

03300 


KK=K+1 

03400 


CK(KK)=0.0 

03500 

03600 


QK(KK)=0.0- 

03700 


-DO 25 IR=2,H 

03800 


R=IR-1 

03900 


CK(KK)=CK(KK)+2.=KAR(IR)=i-C0S(Q*R)) 

04000 


QK<KK)=QK(KK)+BR(IR)=fSIN(Q*R) 

04050 


CPK(KK)=CK(KK) 

04080 


QPK(KK)=OK(KK) 

04100 

04200 

25 

CONTINUE 

04300 


CK(KK)=2.*H*(AR(1 )+CK(KK) + (-1 .=f=*K)*AR(M+1) } 

04400 

04500 


QK(KK)=4.*H+QK(KK) 

04600 

04700 

15 

CONTINUE 

04800 


HH=H+1 

04890- 

04900 


»RITE(33,=f=)(CKa),QKa),I=1,HH) 

05000 


DO 30 1=1 ,HH 
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05100 

05200 

05300 

05400 

05500 

05(400 

05625 

05650 

05700 

05800 

05900 

06000 

06100 

06200 

06300 

06400 

06500 

06600 

06700 

06800 

06900 

07000 

07100 

07200 

07300 

07400 

07500 

07600 

07700 

07B00 

07900 

07950 

08000 

08200 

08300 


QKCX(I) = (0,-1.) 

CKCX<I>=CK(I) 

QKCX(I)=QKCX(I)*QKa) 

GXY(I)=CKCX(I)+QKCX(I) 

30 CONTINUE 

URITE(55,*HCPK(KK) ,QFK(KK) ,KK,KK-1 ,i1M) 
URITE(37,*)(CKCXm,QKCXm,I=1,M}1) 

C SHOOTHING BY HANNING HETHOB 

CK(1)=1CK(1)+CK(2))/2. 
aK(1)=(0K{1)+QK(2))/2. 

BO 37 J=2,i1 

CK(J+1)=(CK(J)+CK(J+1))/2. 

QK(J4-1) = (QK(J)+0K{J+1))/2. 

37 CONTINUE 

DO 40 1=2, H 

CK(I)=(CKa-1)+2.*CK-{I)+CK{I+1))/4. 
QK(I)=(QKa-1)+2.+QK(I)+QKa+1 >)/4. 

40 CONTINUE 

C SHOOTH VERSION 

DO 50 1=1, HN 

CKCX(I)=CK(I) 

QKCX(I)=(0,-1) 

0KCX(I)=OKCX(I)=i=QK(I) 

GXYSa)=CKCX{I)+GKCX(I) 

50 CONTINUE 

■ HRITETJBV* ) ( CKCX ( I ) , QKCXT I ) ; 1=1', )' 

RETURN 

END 



TYPE 
00100 
00200 
00300 
00400 
00500 
00600 
00700 
00800 
00900 
01 000 
01200 
01300 
01400 
01500 
01600 
01700 
01800 
01900 
02000 
02100 
02200 
02300 
02400 
02500 
02600 
02700 
02800 
02900 
03000 
03100 
03200 
03300 
03400 
03500 
03600 
03700 
03800 
03900 
04000 
04100 
04200 
04300 
“0‘4400 
04500 
04600 
04700 
04800 
04900 
05000 
05100 
05200 
05300 
05400 
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PRQGRAH 0IB222 
COHHON BXY(100),6XYS{100) 

COMPLEX FILEX,HT,HS,GXY,6XYS 

DIHENSIOM X(999),Y(999),6KX(100),GKXS(100),GKY{100), 

1 GKYS(100),HT(100),HS(100),C(100),CS(100) 

1 ,D(100),RX(100) 

1 ,A(100),BBnOO),CC(100),C&(100),CSS(100),BC(100) 

REAB(60,*)H 

C OPENS READS, AND CLOSES THE FILE CONTAINING THE 

C DATA FOR THE INPUT FUNCTION, X 

READ(23,=»)HX,NX,TR,(X(I),I=1,NX) 

C CHECKS TO HAKE SURE H, AND N ARE POSITIVE AND EVEN 

C CALCULATES CUTOFF FREQ BANDUIDTH, AND STANDARD ERROR 

CALL CHECK{NX,HX,H,F,B,E) 

URITE(46,50)NX,HX,H,F,B.E 
50 F0RHAT(5X,''NX=',I3,3X,''H='',F15.5,3X,'’H=" 

1 I3,3X,^FC=",F15.5,3X,'B£=%F15.5,3X,^E=%F15.5) 

C READ DATA FOR OUTRPUT FUNCTION, Y 

READ(24,=f)HY,NY,TR,(Y(I),I=1,NY) 

C CHECK y DATA 

CALL CHECK(NY,HY,H,F,B,E) 

C HAKE SURE THERE ARE SAHE NO. OF X AND Y VALUES ■ 

IF<NX .NE. NY)STOP 'NX NOT EQUAL TO NY' 

C AVERAGE X AND Y HADE EQUAL TO ZERO 

C RESULTS OF AVERAGING IN FILE FQR20.DAT 

CALL FORH(X,NX) 

CALL FORH(Y,NY) 

URITE(20,29) 

29 F0RHAT(3X,'X(I)',5X,'Y(I)',5X,'I') 

HRITE(20,31)(X(I),Y(I),I=1,NX) 

31 FORMAT (5X, 13,1 OX, 2E1 0.5, FI 5. 5) 

C POUER SPECTRA FOUND FOR X AND Y 

C RESULTS FOR X FUNCT. F0R30.DAT 

C RESULTS FOR Y FUNCT. F0R31.DAT 

CALL POUER(X,HX,HX,H,BKX,GKXS,n,RX) 

UNITE (30, i-)(GKX(I),GKXS(I), 1=1, M) 

CALL POWER(Y,HY,NY,H,GKY,GKYS,D,RX) 
URITE(31,*)(GKY(I),GKYS(I),I=1,H) 

C COMPLEX VALUED CROSS-SPECTRAL DENSITY FOUND 

C RESULTS IN FILE F0R32.DAT 

CALL CROSS(X,Y,HX,H,NX,A,BB,CC,CD,CSS,DC) 
URITE(32,*)(GXYCI),GXYSa),I=1,M) 

C ''TRANSFER' FONCTlON TO'liND COMPlEX-VALUED 

C COHERENCE FUCTION FOUND. REAL 

CALL TRANS(GKX,GKY,GKXS,GKYS,H,HT,HS,C,CS, 

1 A,BB,CC,CD) 

C FILE IS NAMED FOR STORAGE OF THE RESULTS 

C FILENAME MUST BE 9 CHARACTERS WITH .DAT EXT 

HRITE(46,101) 

101 F0RMAT{12X,'HT',22X,'HS(I)',18X,'C(I)',12X,'CS(I)'///) 
URITE(46,102)(HT(I),HS(I),C(I),CSa),I = 1,M) 

102 FORHAT(2X,2F10.6,3X,2F10.6,3X,F10.6,5X,F10.6///) 

END 
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Computer Program For The Measurement And Analysis Of Shock 
Load And Acceleration P_esPonse(l ) 


Theory 


Dsf initions 


f 

fo 

X & y 


x(t-) 

y(t) 

h 

k 


r 

m 

n 



Tp 

or x^- 

bI 


Ck 


Gk 


'"■k 

G^it) 
^ t 
Gy(f) 

Q 

a5s 

c 


A 

r\j 


Frequency 

l/2h = The cut off frequency- 

variables - Impulse & response wave values 

Record 

Record 

•The time interval between data values 

The harmonic number 

The lag number 

The maximum lag number 

The number of data samples (number of intervals) 

The normalized standard error desired for spectral 
calculations 

The finite analogues of Fourier cosine coefficients 
The finite analogues of Fourier cosine coefficients 
The estimate of the true (autocorrelation function) 
value, R at logr corresponding to the displacement rh 
The period t 
The data values 

The sample standard deviation for x(t) 

The sample stan>dard deviation for y(t) 

The desired equivalent resolution band width- for 
power spectra calculations 

Digital value of the co-'Spectrum estimate at 
harmonic k between x(t) & y(t) 

The "Raw" estimate of the power spectral density 
function at harmonic .k, corresponding to the 
frequency f=kfjj/ra. 

Digital value of the quadspectrum estimate at 
harmonic k between x(t) & y(t) 

Povrer spectral density function of x(t) 

Raw value of power spectrum G;j{-(r) 

Time intervals 

Power"^ spectral density function of y(t) 

Record length 
Index 

Is Parzen Lag weighting function 
Absolute value 

Coherence function at each of the harmonic 
frequency k (f=kfQ/m) 

For smooth 
For raw 
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*^xy ^ ^ ^ 


»y 

Gk 

Gkxs ^ 
GxyTf ) 

G^ys(f ) 

jcvt = E 


X 


The co-spectral density function 

The quadrature spectral density function 

Digital value of the power spectrum estimate at 
harmonic k for x(t) 

Digital value of the power spectrum estimate at 
harmonic k for y(t) 

(G]c»Qk^ ~ G^yCkf^/m) 

The smooth version of G^ 

Gross spectral density function 
Smooth version using Hanning method 
or Hv Time interval for x or y 
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Plc;ltl 2 lru^ Of Contln'uous Pata 

i.e., convertin'? continuous data into discrete numbers. 

The procedure is first to define the points at X'fhich the 
data are observed (sampling). 

If the time interval = At (secs.) between samples = h 
•' sampling rate = 1/h samples per second. 

The useful data will be from 0 to l/2h cps 

The cutoff frequency f^ = l/2h (1) 

knovm as the Nyquist frequency 

For any frequency f in the range 0 ^ f f^, the higher 
frequencies which are aliased with f arej- 

(2fc±f), (4fc+f ), ,(2nfc±f) (2) 

(cos 2Tft = cos 2T{2nfc±f }.l/2fQ = cos iTf/f q) 

Thus all data at frequencies 2nf(,+f have the same cosine 
function as data at frequency f when data are sampled at 
points l/2f^ apart. 

To Deal V7ith Aliasing Problem 

Choose h sufficiently small so that it is nhysically 
unreasonable for data to exist above the associated cutoff 
frequency f^ (example -assume that information below 500 
cps is desired, as the case we are dealing with. This 
means that l/2h -500 or h=l msec, would technically be 
sufficient), ^ 

In general, it Is a good rule to select f«^ Ir oi* 2 
times greater than the maximum frequency of Interest,. 

The second part of concern is quantization, i.e, , the 
actual conversion of the observed values to numerical form. 

If one assumes that the quantization errors follow a 
uniform probability distribution over one scale unit, then 
these errors will have a mean value p(x)=zero and a 
standard deviation of ^^0,29 scale unit. This is the rms 
value of the quantization error, 

For most practical problems,, quantization errors are 
negligible. 

Basic Statistical Calculations For A Single Record 
Digitizing Of Data And Selection Of Sample Size 

Let j n=l,2, ,N be the data values of a record 

U(t) found at the points tp^t^+nh- 

These points are at distance h apart 


J". Ujj=U(tQ+nh) 


( 3 ) 
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- 1 ^ ' 

Sample mean = U = - > Uy, 

W n=l" 


W 


Sample standard deviation = 
And 

For independent data saraules 
= 1/ [/ N 


r 


S = 


N 






( 5 ) 


• The normalized standard error 


Transformation Of Hat a To Zero Mean Value 

This Is done in order that subsequent formulas and 
calculations may be simplified. 

Define a new time history record 
x(t) = U(t) - U , 

Then x(t) has data values ^x^^v given by 

x^ = x{tQ+nh) = - U (6) 

note that x = 0 

The reason for representing the original data values by 
instead of | is to have the S Xj^^ notation indicate 
a zero sample mean -Value, ^ 


Fourier Series Representation 

If a transformed sample record x(t) is periodic with 
a period of Tp and a fundamental frequency f^=l/Tp, then 
x(t) can be represented by the Fourier seriesi 

CO 

x(t) = a^/2 -f H(aq.cos 2irqfit + b^.sin 2‘TTqf^t) 

Where i -- 

a^ = 2/T “ Jx(t) .cos 2'TTqfj^t,dt q= 0,1,2, 

b = 2/T* 1 x(t),sin 2Xqf^t.dt q= 1,2,3,- — 

Assume a sample record x(t) is of finite length 
the fundamental period of the data. 

Also assume that the record Is sampled at an even number 
of N equally spaced points a distance h apart, where h has 
been selected to produce a sufficiently high frequency 
cutoff fQ=l/2h. 

Consider the initial point of the record to be zero and 
denote the transformed data values as before, by 


= x(nh) 


( 7 ) 
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One may calculate the finite version of a Fourier series 
which will pass throu/^h these KT data values. For any point 
t in the interval (0,T^) one obtains t 
A//2. P 

x(t) Aq + Z[A„cos(2Tqt/T ) + II BQain(2"irqt/T ) — (8) 

^ ^ P ^'r! P 

At the particular points t=nh; n=l,2, ,M 


where Tp = Mh 


A'/z 


(V/2)-l 


X = x(nh) = Api + £.Aneos(2T'qn/M) + 21 B_sin(2 irqn/M) 
The coefficients A- & 3 are given byj 

V 4 4 

ll ^ l^n “ ^ ^ ^ 

Aq = -'ll 3C^cos(2 TTqn/N) q=l,2, ,{N/2)-l 

H Y'/ 

“^N/2 ^ ^os TT^n 

= - £2-3ln(2-rqn)/S 

q 1? ^ ^ transformed data,i/herei ^ 

is made equal to zero 


(9) 


( 10 ) 


Autocorrelation And Power Soectra Calculations 

. Basic formulas to calculate autocorrelation and uower 
spectral density functions for single records from 
digitized data 

Sampling Interval 

h = such that h=l/2f ( 11 ) 

where l/f^ is smallest '-'period"*^ in record 

For accurate correlation function measurements i ’.-niere 
the correlation function has frequencies near f., one should 
choose 

If the power spectra measurements are the urlme consider- 
ation, one should choose h=2/5f^. 

Number of . correlation lag values m=l/Be.h (12) 

~Sx Bg=i/ni.h d3) 

Thus Bq will be small for a given h when ra is large 

The sample size ; N is chosen such that N=m/g, ^ (1^) 

The associated minimum record length T^=N,h (-15) 

The number of degrees of freedom for spectral calculations is 
n = = 2M/m (l6) 
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The Normalized Standard Hrror is giTen by £ = I ^ m/N 


^ will be small for a given N when ra is small 


(17) 


Autocorrelation Functions 

r^r If data values a transformed record x(t) 

which is stationary with 2=0, the estimated autocorrelation 
function at the displacement rh is defined by 


A 


A 


Rj, = Hx(rh) = 


/7-r 

H 


n ■ ^n+r 


( 18 ) 


W - r , 1 ^/ 

The Maximum Lag Number m determines the later 
equivalent frequency band -width resolution for the uower 
spectral density function in the frequency interval (0,fc). 
This equivalent band width ist 


Be = 2f^/m = l/(m.h) = l/T'^ax. 

It means B Is twice the range found by dividing the 
frequency interval (0,fj,) into m equally spaced parts 
(fp/m apart). 

(Bo divides the theoretical frequency Interval (-f^^.f^) 
into m equally spaced parts). Thus from knowledge of f^, 
one can choose m in advance so as to have a desired B^. 

Since the quantity m represents the maximum number of 
correlation lag values, the maximum displacement isi 

=">h (20) 

It is desirable to keep the maximum lag ra less than l/lO 
the sample size N. And it is required to smooth the "Raw" 
power spectral density function estimates. 


The autocorrelation function may take on negative as well 
as positive values, 

A normalized value for the autocorrelation function is 
obtained by dividing by R^^ where 

;(0) = i r (21) 


A 

Rr 


A 
= R 


X’ 


A n I A 

I'Hien R^ is normalized, one obtains the quantity R^/Rq 
which theoretically will be between plus and minus one 


-1 1 


( 22 ) 


- S' 
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Power Spectral Hensltj Functions 

For samples data from a transformed record x(t) which 
is stationary with 5=0, a "Raw" estimate ^ true 

power spectral density function is defined for an 

arbitrary f in the ^ranfre 0 f:^ byi 

^^(f) = 2h pQ+2 £^H^.cos(Trrf/f^) + R^. cos { TTnif/f^ )] ( 23 ) 

The total mean square value of the record in the frequency 
ran^e 0 ^ f 


J 


f^ is; 

/V 


G^(f).df = Ro = 


(24) 


It is recommeded that the values of the function G^(f) 
be calculated only at the m+1 special discrete frequencies. 


Where i 


f=k.fp/m k=0,l,2, ,m 


(25) 


2f^/m apai't 


This will provide m/2 Independent spectral estimates 
since speoti’al estimates at points less than 
will be correlated. 

At these discrete frequency polntsi 

/S/ 'v/ 






Smoothing is necessary since the periodo.^>;ram (Raw estimate) 
given by equation (l8) is an inefficient estimate of the 
true spectral density. The variability of these estimates 
does not decrease with Increased record length or sample 
size. 

This leads to the requirement of smoothing the period- 

ogram or equivalently, weighting the correlation function 
nonun iforraly, 

V/hen the computation of frequency response and the 
associated' coherence functions is the main goal, certain 
variations in the above formulas Is appropriate 

■Pp = 1 - 6(r/m)^ + 6(r/m)^ 

= 2(1 - (r/m))^ 

= 0 . 


r=0,l,2, 


•,m/2 


r=(m/2)+l , ,m 

T';^m 


(27) 


Note that Dq = 1 and = 0 

The use of this Parzen Lag 'V;p maintains the sample 
coherence function between its theoretical values of +1, 
as opposed to the Hanning window of equation (28) which 
allows a sample coherence function to vary over wider 
limits depending on the form of the power spectrum. 


The final "Smooth" formula of the power spectral density 
estimate at harmonic k is 

= G^(k.f^/m) = 2h(RQ+ XI n^.Rr* cosClTr.k/m) ) 


where is given by eauatlon 


V = \ 
{'.'■I ^ 
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Cross Correlation Functions 

Choose a maximum las? number ra 
which will give a desired equivalent frequency resolution 
Bg = 1/rah as well as a desired number of decrees of freedom 
2W/m. 

The autocorrelation functions and cower scectral density 
functions are calculated for the transformed variables 
x(t) and y(t) in terms of the data values, separately 

accord inc: to the formulas listed previously. 

The following are the formulas to calculate their joint 
cross-correlation functions and cross-spectral functions 
in terras of the trasformed dataj 


Estimates for the sample cross-correlation functions at 
lag numbers r=0,l,2, ,m are defined by j 


Hxy(rh) 


^^(rh) 


, V-r- 

\ l-i-t 
~ N - r 

n = I 


(29) 

(30) 


(If N is>7m, it may be more convenient to divide by N 
instead of N-r in equations (30) , (29) and (l8),) 


For later determination of the cross-spectral density 
function estimate 
A 

G (f) betvreen x(t) and y(t), calculate for 
r=0,l,2, — — ,m, the two quantities (even and odd parts of 


the cross-correlation function), 

Aj. = A^y(rh) = ^yx^^^O 

^r “ ~ 2 ^^2) 

Cross Spectral Density Functions 

It is a complex-valued quantity defined bys 

*^xy^^) “ ^xy^^^ ” ^ ‘ (33) 


"Raw" estimates from sampled data for the co-spectral 
density function and the quadrature spectral density function 
may be found as follows* 

The formulas are for one-sided spectra which are nonZero 
only for f ^ 0 

At an arbitrary value of f in the range O^f^f^, 

"Raw" estimates for ^ ^ 
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-A '’'"V A 1 

C^y(f) = 2h |_Aq+2 ^^A 3 ,.cos( // rf/f^)+At^,cos(Trmf/fQ^ ^ 34 ^ 

^^y(f) = 2h [^2 sin( irrf/fj,)+B^.sln(ir'Taf'/fQ2] ^ (35) 

Where Siven by equations ( 31 ) ft (32), 

As before, these fimctions should be calculated only at the 
m +1 special discrete frequencies of harmonic. 


At these discrete frequency points* 

= 'c^.yOcf^/ni) =.2h [A^+2t AyCos(TTrk/m)+(-^ ^3^^ 

I I “ I 

% = ' 

^ \'=\ 

’'Smooth estimates for at harmonic k may hou be 

calculated as in equation ( 2 ? ) , by usinf; the Hannin.*? method. 

This yield* 



0.5Cq + 0.5Cj 
0.5Qo + 0.5Qi 
0.25Cj^._^ + 0.5C]^ + 
+ 0.5Qj, + 
o.icji.i + o.5c„ 
0.5%_i + o.5Q„ 


0.25C]J+1 

0 . 254+1 



m) 


At the m+1 discrete frequencies one obtains the smooth" 
cross-spectral density estimates 

Frequency Response Functions 


H(f) is a complex-valued function 

H(f ) = I H(f ) I ,e-J25(f) ■ 

• '^-•^bere * - 

(H{f)) is the gain factor of the system 
0(f) is the phase factor of the system 
If a stationary input x(t) produces an output y(t), 
the system gain factor 

|H(f)| = (Gy{f)/Cr3 +(f))‘- 


(39) 


(^^ 0 ) 
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Hence, at the discrete frequencies f=kfc/ra, the ^aln 
factor may be estimated digitally by 

l^kl = )H(kf^/m)) = 

For either single input -single output problems where, 
extraneous noise is present only at the output, or multiple 
input— single output problems vrhere the inuuts are 
uncorrelated, another more general method' for estimating 
the system frequency response function (including both gain 
and phase factors) is given byj- 

H(f) = J^y(f)/S’^(f) — ..(^1-2) 

It follows thatj- 

|H(f)| = |S;jy(f)| A^lt) 

And A A 

0(f) = Q^y(f) 

Therefore at the discrete frequencies f='kf^/m the gain 
factor and phase factor may be estimated digitally byj- 
, A • A9 Ap 1 A 

|Hkl = ^ (44) 

^ ^ A A ' 

0 ^ = tan‘*^(Q^/Cj^) (45) 



Coherence Functions 


The coherence function 


2 

between two stationary 


records x(t) and y(t) is defined byi- 


if) 

xy 


1g (f ) \ /(G^(f).G (f)) 


'zy' ' ' " y' 

The coherence function theoretically should satisfy 


(46) 


0 ^ ^ ^ 


for all f. 


In terras of digital calculations, at the discrete 
frequencies f=lcf^/m the coherence function is estimated by 
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The Transfer function is the inatheraatical description 
of the system. It can also be used to measure the relat- 
ionship between any two sisals. It can be defined asj- 

Average cross power spectrum of input ^output 

Transfer Function = — 

Average power spectrum of input 
This gives more reliable transfer function. 

The coherence function can be used to check the validity 
of the transfer function, i. e. , it can be used, to measure 
the degree of causality between any two signals. A "0" 
value means no coherence between input and output of the 
system. A value of "1” means there is a complete coherence 
between them indicating that there is only one input and 
the system is -linear. 



